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Abstract 



We discuss the most important Physics thus far extracted from studies of B meson decays. 
Measurements of the four CP violating angles accessible in B decay are reviewed as well as direct 
CP violation. A detailed discussion of the measurements of the CKM elements Vet and Vub from 
semileptonic decays is given, and the differences between resulting values using inclusive decays 
versus exclusive decays is discussed. Measurements of "rare" decays are also reviewed. We point 
out where CP violating and rare decays could lead to observations of physics beyond that of the 
Standard Model in future experiments. If such physics is found by directly observation of new 
particles, e.g. in LHC experiments, B decays can play a decisive role in interpreting the nature of 
these particles. 
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1 Introduction 

The forces of nature generally reveal their properties by how they act on matter. On the most 
fundamental material scale that we are aware of matter is formed from fermions. These take 
two forms, leptonic matter and quark matter. The former do not have any strong interactions, 
that is they lack a property called "color charge" , which allow quarks to bind together either 
mesonically or baryonically. New and therefore as yet unknown forces could effect both leptons 
and quarks. Here, we concentrate on how such forces effect quarks, especially the b quark. 

Light matter consists mostly of u, d and s quarks. In 1963 Cabibbo showed that weak 
interactions of mesons and baryons containing s quarks were suppressed with respect to those 
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without s quarks by an amount tan^c", where the "Cabibbo" angle 9c must be determined 
experimentally [1]. The s was further shown to have an important and at that time a 
mystifying role, by the discovery of CP violation in decays in 1964 [2,3]. When the c quark 
was discovered in Nov. 1974 [4] (though its existence was speculated earlier [5]), it became clear 
that 6c was the mixing angle between two quark doublets, (c, s) and {u, d). However, it is 
not possible to generate a CP violating phase with only two quark doublets. 

This was recognized even before the discovery of the charm quark by Kobayashi and Maskawa, 
who postulated the existence of yet another quark doublet (6, t) [6], in work for which they 
were awarded the Nobel Prize in 2008. While the t quark is the heaviest, having a mass of 173 
GeV, they are difficult to produce and decay before they can form a hadron, thus excluding 
many otherwise possible studies. Much interesting work has been done with the s and c 
quarks, but in this article we discuss the physics of the b quark, which turns out to be the 
most interesting of the six quarks to study. 

1.1 How B's Fit Into the Standard Model 

First we will discuss how particles formed with 6-quarks fit into current paradigm of particle 
physics, the "Standard Model" (SM) [8]. The SM has at its basis the gauge group 
SU(3)xSU(2)xU(l). The first term corresponds to the strong interaction and SU(3) describes 
the octet of colored gluons which are the strong force carriers of quantum chromodynamics. 
SU(2)xU(l) describes the weak interaction and is the product of weak isospin and 
hypercharge. We speak of the fundamental objects being spin-1/2 quarks and leptons and the 
force carriers generally spin-1 objects. The spin-0 Higgs boson, yet to be discovered, is 
necessary for generating mass [9]. 

Particles containing h quarks can be 5°, Bs, or B^ mesons, depending on whether the 
light anti-quark that it pairs with is d, u, s, or c, or a baryon containing two other quarks. 
Mesons containing b and b quarks are also interesting especially for studies of quantum 
chromodynamics, but will not be discussed further in this article, as we will concentrate on 
weak decays and discuss strong interactions as an important and necessary complication that 
must be understood in many cases to extract information on fundamental b quark couplings. 

The quarks come in three repetitions called generations, as do the leptons. The first generation 
is d u, the second s c and the third b t. In the second and third generations the charge +2/3 
quark is heavier than the charge -1/3; the first generation has two very light quarks on the 
order of a few MeV with the d thought to be a bit heavier^ Decays usually proceed within 

^Isospin invariance is related to the equality of it and d quark masses. The PDG [7] gives the u quark mass 
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generations, so the c decays predominantly to the s quark via the quark level process 

c — > W^s, though some decays do go to the first generation as c — > W^d. The ratio of these 

amplitudes approximate the Cabibbo angle discussed earlier. 

The mixing matrix proposed by Kobayshi and Maskawa [6] parameterizes the mixing between 
the mass eigenstates and weak eigenstates as couplings between the charge +2/3 and -1/3 
quarks. We use here the Wolfenstein approximation [10] good to order A^: 

/ Vud 



VcKM 




(1) 



A A\^{p - IT]) 

1 - AV2 - A^(l + 4^2)/8 
-AX^ + A\\l/2-{p + ir])) l-A^\^/2 

In the Standard Model A^ A, p and ry are fundamental constants of nature like G, or uem] i] 
multiplies the imaginary i and is responsible for all Standard Model CP violation. We know 
A=0.226, A ~0.8 and we have constraints on p and r]. Often the variables p and rj are used 
where 

p + ifi={p + ir,){l - AV2) = (2) 
where the definition in terms of the CKM matrix elements is correct to all orders in A [11]. 

Applying unitarity constraints allows us to construct the six independent triangles shown in 
Figure [H Another basis for the CKM matrix are four angles labeled as x' two of a, 

f3 and 7 since a + /? + 7 = vr [12]. (These angles are also shown in Figure [TJ) CP violation 
measurements make use of these ang lesi 

B meson decays can occur through various processes. Some decay diagrams with intermediate 
charged vector bosons are shown in Figure [21 The simple spectator diagram shown Figure [2]^a) 
has by far the largest rate. Semileptonic decays proceed through this diagram, and allow us to 
measure the CKM couplings Vcb and Vub by considering only the hadronic uncertainties due to 
the spectator quark. The color suppressed diagram Figure WIJo) exists only for hadronic decays. 
It can occur only when the colors of the quarks from the virtual W~ decay match those of the 
initial B meson. Since this happens only 1/3 of the time in amplitude, the rate is down by 
almost an order of magnitude from the spectator decays. The annihilation Figure [2]^c) 
describes the important decay B^ t~V and will be discussed in detail later. The W 
exchange Figure [2||d) diagram is small. The box diagram Figure [2^e) is the source of mixing in 



between 1.5-3.3 MeV and the d mass between 3.5-6.0 MeV, where the large range indicates the considerable 
uncertainties. 

^The Belle collaboration defines cf>2 = a, 4>i = P, and (^3 = 7. 
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Figure 1: The 6 CKM triangles resulting from applying unitarity constraints to the indicated 
row and column. The CP violating angles are also shown. 
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Figure 2: Some B decay diagrams. 



the B system. An analogous diagram exists also for the Bg meson. Finally the Penguin 
diagram Figure [2)[f) is an example of one of several loop diagrams leading to "rare decays" 
that will be discussed in detail in a subsequent section. 



1.1.1 Dark Matter 

"Dark Matter" was first shown to exist by Zwicky studying rotation curves of galaxies [13] . 
The motion could only be explained if there was massive cloud of matter that was not 
luminous. We still do not know what composes this dark matter, though hopes are it will be 
discovered at the LHC. An even more mysterious phenomena called "Dark Energy" may also 
have a connection to particle physics experiments [14], perhaps via "Extra Dimensions" [15]. 



6 



1.1.2 Baryogenesis 

When the Universe began with the Big Bang, there was an equal amount of matter and 
antimatter. Now we have mostly matter. How did it happen? Sakharov gave three necessary 
conditions: Baryon (B) number violation, departure from thermal equilibrium, and C and CP 
violation [16]. (The operation of Charge Conjugation (C) takes particle to anti-particle and 
Parity (P) takes a vector to — T*.) 

These criteria are all satisfied by the Standard Model. B is violated in Electroweak theory at 
high temperature, though baryon minus lepton number is conserved; in addition we need 
quantum tunneling, which is powerfully suppressed at the low temperatures that we now have. 
Non-thermal equilibrium is provided by the electroweak phase transition. C and CP are 
violated by weak interactions. However the violation is too small. The ratio of the number of 
baryons to the entropy in the observed part of the Universe needs to be ~ 5 x 10^^^, while the 
SM provides many orders of magnitude less. Therefore, there must be new physics [17]. 

1.1.3 The Hierarchy Problem 

Our worry is why the Planck scale at ~ 10^^ GeV is so much higher than the scale at which we 
expect to find the Higgs Boson, ~100 GeV. As Lisa Randall said [18] "The gist of it is that the 
universe seems to have two entirely different mass scales, and we don't understand why they 
are so different. There's what's called the Planck scale, which is associated with gravitational 
interactions. It's a huge mass scale, but because gravitational forces are proportional to one 
over the mass squared, that means gravity is a very weak interaction. In units of GeV, which 
is how we measure masses, the Planck scale is 10^^ GeV. Then there's the electroweak scale, 
which sets the masses for the W and Z bosons. These are particles that are similar to the 
photons of electromagnetism and which we have observed and studied well. They have a mass 
of about 100 GeV. So the hierarchy problem, in its simplest manifestation, is how can you have 
these particles be so light when the other scale is so big." We expect the explanation lies in 
physics beyond the Standard Model [19]. 

1.2 B Decays as Probes for New Physics 

When we make measurements on B decays we observe the contributions of SM processes as 
well as any other processes that may be due to physics beyond the SM or New Physics (NP). 
Other diagrams would appear adding to those in Figure [2] with new intermediate particles. 
Thus, when it is declared by those who say that there isn't any evidence of NP in B decays, we 
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have to be very careful that we have not absorbed such new evidence into what we declare to 
be SM physics. There are several approaches that can be followed. 

One approach is to simply predict the decay rate of a single process in the SM with known 
couplings and compare to the measurements. The classical case here is 6 — > S7 and we will 
discuss this and other specific examples later. Another approach is make different 
measurements of the CKM parameters in different ways and see if they agree. This is normally 
done by measuring both angles and sides of the CKM triangle, but other quantities can also be 
used. This is the approach used by the CKM fitter [11] and UT fit groups [20]. In yet a third 
approach, the exact same quantity can be measured in several ways, even if cannot be 
predicted in the SM. An example here is measuring the CP violating angle (3 using 

— > J/ijjKs decays that proceed through the diagram in Figure [2][b), at least in the SM, and 
another process that uses the "Penguin" diagram in Figure [2)||f), e.g. — > 4>Ks- 

The punch line is that if new, more massive particles exist in a mass range accessible to the 
LHC then they MUST contribute to rare and CP violating B decays! Even if measurements 
are precise enough only to limit the size of these effects, the properties of these new particles 
will be much better understood. This is the raison d'etre for the further study of B decays. 

2 Measurements of Mixing & CP Violation 
2.1 Neutral B Meson Mixing 

Neutral B mesons can transform into their anti-particles before they decay. The diagrams for 
this process are shown in Figure [3] for the B^- There is a similar diagram for the Bg. Although 
u, c and t quark exchanges are all shown, the t quark plays a dominant role mainly due to its 
mass, as the amplitude of this process is proportional to the mass of the exchanged fermion. 



Under the weak interactions the eigenstates of flavor degenerate in pure QCD can mix. Let the 
quantum mechanical basis vectors be {|1), |2)} = {l-B*^), \B^)^. The Hamiltonian is then 




Figure 3: The two diagrams for B^ mixing. 
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2 I M 

The Schrodinger equation is 




B\t)) I \B^(t)) 




n = M-^T= )-,( - ). (3) 



B'{t)) \ ^ ( I B^{t)) 



W (4) 



Diagonalizing we have 

Am = — ms^ = 2 IM12I (5) 

Ar = Tl-Th = 2\Ti2\cos^ (6) 

where H refers to the heavier and L the hghter of the two weak eigenstates, and 

(f) = arg (— Mi2/ri2). We expect that AF is very small for mesons but should be significant 

for Eg mesons. 

Bd mixing was first discovered by the ARGUS experiment [21] (There was a previous 
measurement by UAl indicating mixing for a mixture of and B^ [22]). At the time it was 
quite a surprise, since the top-quark mass, m^, was thought to be in the 30 GeV range. Since 
6-flavored hadrons are produced in pairs, it is possible to observe a mixed event by having both 
-B's decay semileptonically. Thus the number of events where both Bs decay into leptons of 
the same sign is indicative of mixing. We can define R as the ratio of events with same-sign 
leptons to the sum of same-sign plus opposite-sign dilepton events. This is related to the 
mixing probability. The OPAL data for R are shown in Figure [4] [23]. 

Data from many experiments has been combined by "The Heavy Flavor Averaging Group," 
(HFAG) to obtain an average value Amd = (0.507 ± 0.004) x 10^^ pg-i [24]. 

The probability of mixing is related to the CKM matrix elements as [3, 25] 

_ Am G]p 2 |T/*T/ |2 2tp( 

Xd = — = -^BBjB'mBTB\Vt^ytd\ ^tFyjfT jVQCD, (7) 

where Bb a parameter related to the probability of the d and b quarks forming a hadron and 
must be estimated theoretically, F is a known function which increases approximately as m^, 
and riQCD is a QCD correction, with value about 0.8. By far the largest uncertainty arises 
from the decay constant, /b- 

In principle /b can be measured. The decay rate of the annihilation process B^ i~V is 
proportional to the product of /^l^ufep. Of course there is a substantial uncertainty associated 
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Figure 4: The ratio, R, of same-sign to total dilepton events as a function of proper decay 
time, for selected B — > D*~^Xi~v' events. The jet charge in the opposite hemisphere is used to 
determine the sign correlation. The curve is the result of a fit to the mixing parameter. 



with \Vub\- The experimental evidence for B~ t~V is discussed in Section [STTl and 
substantial uncertainty exists in the branching ratio measurement. Thus, we need to rely on 
theory for a value of fs- 

The ratio of Bg to B^ mixing frequency, however, provides a better situation in terms of 
reducing model dependent errors. Using Eq. [7] for B^ mixing and an analogous relation for Bg 
mixing and then dividing them results in 



Xd ^ Bb fl TUB TB \VtlVtd\'^ 
Xs ~ Bb, fl TUB, TB, \ytb^ts\^ 



(8) 



The CKM terms are 



V:,Vtd\' = A\'\{l-p-i^)\' = AX-'{p-lY + r^' ^nd (9) 



ignoring the higher order term in Vts- Solving the ratio of the two above equations for |^ 
gives a circle centered at (1,0) in the p — rj plane whose radius depends on Xd/xg- The 
theoretical errors are now due to the difference between having a light quark versus a strange 
quark in the B meson, called SU(3) splitting. 

For many years experiments at the using e^e~ colliders at both LEP and the SLC had set 
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lower limits on B'^ mixing [7]. In 2006 Bg mixing was measured by the CDF collaboration [26]. 
The DO collaboration had previously presented both a lower and an upper limit [27]. 

We now discuss the CDF measurement. The probability, 'P(t) for a Bg to oscillate into a Bg is 
given as 

r{t) (S, ^ Bs) = ^r,e-^^* [1 + cos {Amst)] , (10) 
where t is the proper time. 

An amplitude A for each test frequency uj, is defined as [28] 

V{t) = ^Tge-^^' [1 + A cos {ujt)] . (11) 

For each frequency the expected result is either zero for no mixing, or one for mixing. No other 
value is physical, although measurement errors admit other values. Figure [5] shows the CDF 
results. 




5 10 15 20 25 30 

Arrig [ps ^] 

Figure 5: The measured amplitude values and uncertainties versus the BgBg oscillation fre- 
quency Anis for a combination of semileptonic and hadronic decay modes from CDF. 



At Arus = 17.75 ps~^, the observed amplitude A = 1.21 ±0.20 (stat.) is consistent with unity, 
indicating that the data are compatible with BgBg oscillations with that frequency, while the 
amplitude is inconsistent with zero: A/ga = 6.05, where is the statistical uncertainty on A, 
the ratio has negligible systematic uncertainties. (This then is called a "60" effect.) The 
measured value is 

Arus = 17.77 ± 0.010 ± 0.07 ps"^ , (12) 
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where the first error is statistical and the second systematic. 

In order to translate the mixing measurements to constraints on the CKM parameters p and r/, 
we need to use theoretical values for the ratios and in Eq. [HI It is interesting that in 
practice it is not possible to measure either of these numbers directly. It is usually assumed, 
however, that which could in principle be measured via the process B~ — > t~v is the 
same as fs- There is no way to measure fss- The charm system, on the other hand, provides 
us with both i'^u and i'^u, and both rates have been measured. They also have 

been calculated in unquenched lattice quantum electrodynamics (QCD). 

The combined efforts of the HPQCD and UKQCD collaborations predict fD+ = (207 ± 4) 
MeV [29], while the CLEO measurement is in astonishing agreement: (205.8ib8.5ib2. 5) 
MeV [301. Furthermore, the measurements of fr^+ are not in such good agreement with the 
Follana et al. calculation of (241ib3) MeV. The average of CLEO and Belle results as 
determined by Rosner and Stone is (273ibl0) MeV [31]. The discrepancy is at the 3 standard 
deviation level and not yet statistically significant, though it bears watching. 

Unfortunately, the group that has calculated /£>+ is not the same group that has determined 
fB+- The theoretical calculations used for the Bb terms and the Jb terms are summarized by 
Tantalo [32]. He suggests using values of /^^ = (268 ib 17 ib 20) MeV, 

fBjfB = 1.20 lb 0.02 lb 0.05, Bb, = 0.84 ib 0.03 ib 0.055 and Bb = 0.83 ib 0.01 ib 0.06. These 
numbers allow us to measure the length of one side of the CKM triangle using Eq. [H Use of 
this measurement will be discussed in more detail in Section [6l 

2.2 CP Violation in the B System 

We have two quantum mechanical operators: Charge Conjugation, C, and Parity, P. When 
applied to a particle wavefunction C changes particle to antiparticle and vice-versa. Applying 
P to a wavefunction '(/'(r) we have Pip(r) = ■ip{—r). The P operator can be thought of 
changing the natural coordinate system from right-handed to left-handed. If nature was blind 
to handedness, then P would be always conserved. By applying the P operator twice we end 
up with P^'i/;(r) = 'ip{r), so the eigenvalues of P are ±1. Therefore wave-functions, or particles 
represented by such wave-functions, have either intrinsic positive parity +1 (right-handed) or 
-1 (left-handed). 

Weak interactions, characterized by a combination of vector minus axial-vector currents, are 
known to be left-handed. Therefore, handedness matters and its well known that Parity is 
maximally violated in weak decays [3]. Since C changes left-handed particles to right-handed 
anti-particles, the product CP symmetry could have been preserved, but nature decided 
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otherwise. Different particle transitions involve the different CP violating angles shown in 
Figured) Measuring these independently allows comparisons with measurements of the sides of 
the triangle and any differences in constraints in p and rj can be due to the presence of new 
physics. 

Consider the case of a process B ^ f that goes via two amplitudes, A and B each of which has 
a strong part e. g. s_a and a weak part w^. Then we have 

r{B^f) = (|^|e^('^+"'^) + |e|e^('e+"'B))^ (13) 

T(B^J) = (|^|e^(^'4-"''4) + |^|e^(^e-"'e))^ (14) 
r{B ^ f)-T(B ^J) = 2\AB\sm{sA-SB)sm{wA-WB) • (15) 

Any two amplitudes will do, though its better that they be of approximately equal size. Thus 
charged B decays can exhibit CP violation as well as neutral B decays. In some cases, we will 
see that it is possible to guarantee that |sin(syi — ss)| is unity, so we can get information on 
the weak phases. In the case of neutral B decays, mixing serves as the second amplitude. 



2.2.1 Formalism of CP Violation in Neutral B Decays 

Consider the operations of Charge Conjugation, C, and Parity, P: 

C\B{jf)) = \B{lf)), C\B{t)) = \B{t)) (16) 

P\Bit)) = -\Bi-lf)), P\Bit)) = -\Bi-t)) (17) 

CP\B{t)) = -\B{-t)), CP\Bit)) = -\Bi-t)) ■ (18) 
For neutral mesons we can construct the CP eigenstates 

\B',) = -^{\B')-m) , (19) 

= l=[\B^) + \-^)) , (20) 

where 

CPIB?) = , (21) 

CP\B^2) = -\B'2) ■ (22) 

Since B^ and B^ can mix, the mass eigenstates are a superposition of a\B^) + b\B^) which 
obey the Schrodinger equation 
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If CP is not conserved then the eigenvectors, the mass eigenstates \Bi) and are not the 

CP eigenstates but are 

\BL)=p\B^) + q\'^)^ \BH)=p\B^)-q\'^)^ (24) 



where 

1 l + es 1 l-es 



" V2Vi + ki?P' ^/2v/rTk^' 

CP is violated if 7^ 0, which occurs if / 1. 



(25) 



The time dependence of the mass eigenstates is 

\BL{t)) = e-^^*/2gimLt/2|^^(Q^^ (26) 
\BH{t)) = e-r«*/2gimHi/2|5^(o)), (27) 

leading to the time evolution of the flavor eigenstates as 

\B\t)) = e-(^™+^)*fcos^|P°(0))+z^sm^|#(0))) (28) 

V 2 p 2 J 

|P°(t)) = e-(-+i)*fz^sin^|pO(0))+cos^|P°(0))), (29) 

\ q 2 2 / 

where m = {mi + mii)j2^ Am = mH — mi and F = ~ F/^, and t is the decay time in the 
B^ rest frame, the so-called "proper time" . Note that the probability of a B^ decay as a 
function of t is given by {B^ {t)\B'^ {t))* , and is a pure exponential, e"'"*/^, in the absence of CP 
violation. 



2.2.2 CP Violation for B Via the Interference of Mixing and Decay 

Here we choose a final state / which is accessible to both B^ and b'^ decays [3]. The second 
amplitude necessary for interference is provided by mixing. Figure [6] shows the decay into / 
either directly or indirectly via mixing. It is necessary only that / be accessible directly from 
either state; however if / is a CP eigenstate the situation is far simpler. For CP eigenstates 

CP|/cp) = Mfcp). (30) 



It is useful to define the amplitudes 

A = {fcp\n\B"), A={fcp\n\W). (31) 



If 



A 
A 

detai 



7^ 1, then we have "direct" CP violation in the decay amplitude, which we will discuss in 
later. Here CP can be violated by having 

A = ^ • 4 ^ 1, (32) 
p A 
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Figure 6: Two interfering ways for a to decay into a final state /. 

which requires only that A acquire a non-zero phase, i.e. |A| could be unity and CP violation 
can occur. 

Other useful variables, that are independent of any phase convention are 

(pu = arg ' ^^'^ ^^^^ 

T^2 1 - \q/p\' _ , 

The first quantity can be related to CKM angles, while the second can be measured by the 
"semileptonic asymmetry," or for that matter in any flavor specific decay [33]: 

, T[B°(t) ^ rX]-T[B°(t) ^ £+X] AT 

As Lit) = = -rjj tan c^u , 34 

for either B^ or for Bs mesons, separately. 

A comment on neutral B production at e~^e~ colliders is in order. At the T(4S') resonance 
there is coherent production of B^B^ pairs. This puts the B^s in a C = —1 state. In hadron 
colliders, or at e^e~ machines operating at the Z^, the B^s are produced incoherently. For the 
rest of this article we will assume incoherent production except where explicitly noted. 

The asymmetry, in the case of CP eigenstates, is defined as 

^ r{B^{t)^fcp)-r(B'{t)^fcp) 
''^''^ ~ r [B\t) - fcp) + r - jcp) ' 

which for |g/p| = 1 gives 

(1 - |Ap) cos (Ami) - 2ImAsin(Amt) 

«/cp = rT]A|2 ■ ^ -* 

For the cases where there is only one decay amplitude |A| equals 1, and we have 

aj^p = — ImAsin(Ami). (37) 
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Only the factor — ImA contains information about the level of CP violation, the sine term is 
determined by B mixing. In fact, the time integrated asymmetry is given by 

^fcp = ~YT^^^^ " -0.48ImA . (38) 

This is quite lucky for the study of mesons, as the maximum size of the coefficient for any x 
is —0.5, close to the measured value of Xd- 

Let us now find out how ImA relates to the CKM parameters. Recall ^ = p ' The first term 
is the part that comes from mixing: 

= = , (1^4^M! ^ = e-^^^ and (39) 

l^'- = - '^'-/\ =su.m. (40) 
P (l-p)^ + r/2 



2.2.3 Measurements of sm{2(]) 

To evaluate the decay part we need to consider specific final states. For example, consider the 
final state J/ ipKs ■ The decay diagram is shown in Figure [71 In this case we do not get a phase 
from the decay part because 



A \VcbVc 

is real to order 1/A^. 



cs I 



(41) 



In this case the final state is a state of negative CP, i.e. CP\J /tpKs) = —{J/ipRg). This 
introduces an additional minus sign in the result for ImA. Before finishing discussion of this 
final state we need to consider in more detail the presence of the Kg in the final state. Since 
neutral kaons can mix, we pick up another mixing phase (similar diagrams as for B^, see 
Figure [3|) . This term creates a phase given by 



P^K \VcdVc 



cs I 



which is real to order A^. It necessary to include this term, however, since there are other 
formulations of the CKM matrix than Wolfenstein, which have the phase in a different 
location. It is important that the physics predictions not depend on the CKM conventionlfl 



^Here we don't include CP violation in the neutral kaon since it is much smaller than what is expected in the 
B decay. 
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In summary, for the case of / = J/ij^Ks^ ImA = — sin(2/3). The angle (3 is the best measured 
CP violating angle measured in B meson decays. The process used is — > J/ipKs, although 
there is some data used with the ip{2S) or with K]^ (where the phase is +sin(2/3)). 




Figure 7: The Feynman diagram for the decay B^ J/ipK^. 

Although it is normally thought that only one decay amplitude contributes here, in fact one 
can look for the presence of another source of CP violation, presumably in the decay 
amplitude, by not assuming |A| equals one in Eq. [3^1 Then the time dependence of the decay 
rate is given by 

2ImA . , 1 - IAP , 

""fcp = ~TT\W ^'""^ ""^^ " i + |A|2 cos(Amt). (4 

Thus aff-,p{t) has both sin(Amt) and cos(Amt) terms, and the coefficients of these terms can 
be measured. Let us assign the labels 



^~ i+iAr i+iA|2- ^^^^ 

(Note that the sign of the S term changes depending on the CP oi the final state, but not the 
C term.) 

The most precise measurements of sin(2/3) have been made by the BaBar and Belle 
experiments [34]. These measurements are made at the T(4S') resonance using 
e~^e~ T(45') B^B^ , and with one of the neutral i?'s decaying into J/tpK^. Because the 
T(4S') is a definite quantum state with C = — 1, the formulae given above have to modified 
somewhat. One important change is the definition of t. The time used here is the difference 
between the decay time of the J /iIjK^ and the other B^ , also called the "tagging" B because we 
need to determine its flavor, whether B^ or B^ , in order to make the asymmetry measurement. 

While we will not discuss flavor tagging in general, it is an important part oi CP violation 
measurements. At the T(45') once such very useful tag is that of high momentum lepton as a 
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6-quark, part of a -B meson decays semileptonically into an £~ , while a 6-quark decays into an 
e^. Other flavor tagging observables include the charge of kaons and fast pions. Modern 
experiments combine these in a "neural network," to maximize their tagging efficiencies [35]. 
The BaBar data are shown in Figure [8l Belle has similar results. 




-5 5 

At (ps) 



Figure 8: (a) The number of J/ipKs candidates in the signal region with either a tag (N^o), 
or a B tag (N—o) as a function of At. (b) The raw asymmetry, (A''^o — N—o)/{Nqo + -^-go). 

The solid (dashed) curves represent the fit projections as functions of At for both B^ and 5^ 
tags. The shaded regions represent the estimated background contributions. 

The average value as determined by the Heavy Flavor Averaging Group is 

sin(2/3) = 0.671 it 0.024, where the dominant part of the error is statistical. No evidence is 

found for a non-zero C term with the measured average given as 0.005ib0.020 . 

Determining the sine of any angle gives a four-fold ambiguity in the angle. The decay mode 
B^ J/tf^K*^, K*^ — > Kgir^ offers a way of measuring cos 2/3 and resolving the ambiguities. 
This is a subtle analysis that we will not go into detail on [36]. The result is that 
P = 2imtlit degrees. 

2.2.4 Measurements of a 

The next state to discuss is the CP+ eigenstate / = 7r"'"7r~ . The simple spectator decay 
diagram is shown in Figure (Ha). For the moment we will assume that this is the only diagram, 
though the Penguin diagram shown in Figure O^b) could also contribute; its presence can be 
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Figure 9: Tree (a) and Penguin (b) processes for neutral B decay into either vt^tt or p+p . 




-0.8 - i - 

-0.8 -0.6 -0.4 -0.2 

s 

Figure 10: Coefficients of the sine term S and the cosine term C in time dependent CP violation 
for neutral B decay into tt^tt", showing BaBar and Belle results, and the HFAG average. 
Contours are shown at 60.7% confidence level. 



inferred because it would induce a non-zero value for C, the coefficient of the cosine term in 
Eq. [321 For this b — > uud process we have 



and 



e-"^'^ (45) 



Im(A) = Im(e-2*/3e-2n) = ini(e2^°) = - sin(2a) (46) 



Time dependent CP violation measurements have been made by both the BaBar and Belle 
collaborations [37]. Both groups find a non-zero value of both C and 5, though their values are 
not in particularly good agreement. The HFAG average is shown in Figure [TOl along with the 
experimental results. The value of C clearly is not zero, thus demonstrating direct CP 
violation. (Historically, this was an important observation because it showed that CP violation 
could not be due to some kind of "superweak" model ala Wolfenstein [38].) 
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Table 1: Branching ratios B ^ pp modes in units of 10 



Mode 


BaBar 


Belle 


Average [7,24] 


p+p- 




22.8 ± 3.8l^:^ 


24.2 ±3.1 


p+p^ 


16.8 ±2.2 ±2.3 


31.7±7.lt|? 


18 ±4 


p^p^ 


0.92 ±0.32 ±0.14 


0.4l[j:| ± 0.3 


n 74+0-30 

'-'•'^-0.27 



The non-zero value of C shows the presence of at least two amplitudes, presumably tree and 
penguin, in addition to the mixing amplitude, making extraction of a difficult. All is not lost, 
however. Gronau and London [39] showed that a could be extracted without theoretical error 
by doing a full isotopic spin analysis of the tttt final state. The required measurements include: 
C , S, B{B+ 7r+7r°), ^ 7r°7r°), and vrOyrO). The last two items require a 

flavored tagged analysis that has not yet been done and whose prospects are bleak. Grossman 
and Quinn have showed, however, that an upper limit on T(B^ — > 7r'^7r^)/r(i?"^ — > vr+vr^) can 
be used to limit the penguin shift to a [40]. Unfortunately, current data show a relative large 
7r°7r°) = (1.55 ± 0.19) x 10^^ rate, compared with 

= (5.16 ± 0.22) X 10-6, implying a large penguin contribution [7, 24] , and the 
limit is very weak. 

Use of the p^ p^ final state is in principle quite similar to vr+vr^, with some important caveats. 
First of all, it is a vector-vector final state and therefore could have both CP+ and CP- 
components. It is however possible, doing a full angular analysis to measure the CP violating 
phase separately in each of these two amplitudes. The best method for this is in the 
"transversity" basis and will be discussed later [41] in Section 12.2.61 It is possible, however, for 
one polarization component to be dominant and then the angular analysis might not be 
necessary. In fact the longitudinal polarization is dominant in this case. BaBar measures the 
fraction as 0.992 ± 0.024l°;°26 [42], and Belle measures it as 0.94i;°;°|^ ± 0.030 [43]. Thus we 
can treat this final state without worrying about the angular analysis and just apply a 
correction for the amount of transverse amplitude. 

In addition, p mesons are wide, so non-i? backgrounds could be a problem and even if the 
proper B is reconstructed, there are non-resonant and possible oivr contributions. 
Furthermore, it has been pointed out that the large p width could lead to the violation of the 
isospin constraints and this effect should be investigated [44]. The relevant branching ratios 
are given in Tabled) 

Nevertheless, the small branching ratio for p^ p^ , if indeed it has been observed at all, shows 
that the effects of the Penguin diagram on the extracted value of sin(2a) are small, and this 
may indeed be a good way to extract a value of a. The time dependent decay rates separately 
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Figure 11: The number of p candidates in the signal region with either a tag (a) (N^o), 



or a tag (b) (-^-go) as a function of At. (c) The raw asymmetry, (A''^o — N^))/{N^o + A''-^ 
The dashed curves represent the estimated backgrounds. 



for P^P^ and B — > p^p^ and their difference are shown in Figure [TT] from BaBar. In is 

interesting to compare these results with those in Figure El We see that the measured 
asymmetry in the p^p^ decay more or less cancels that from B^ — B^ mixing inferring that 
sin(2Q;) is close to zero. 

Results from the time dependent CP violation analysis from both Belle and BaBar are shown 
in Figure [T2l 

The data can be averaged and a determined by using the isospin analysis and the rates listed 
in Table [H Unfortunately the precision of the data leads only to constraints. These have been 
determined by both the CKM fitter group [45] and the UT fit group [20]. These groups 
disagree in some cases. The CKM fitter group use a frequentist statistical approach, while UT 
fit uses a Bayseian approach. The basic difference is the that the Bayesian approach the 
theoretical errors are taken as having a Gaussian distribution. Here we show in Figure [13] the 
results from CKM fitter. 

The final state vr+vr^Tr'^ can also be used to extract a. Snyder and Quinn proposed that a 
Dalitz plot analysis of — > pvr ^ vr^vr^vr'' can be used to unravel both a and the relative 
penguin-tree phases [46]. The Dalitz plot for simulated events unaffected by detector 
acceptance is shown in Figure [HI 
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Figure 12: Coefficients of tlie sine term S and the cosine term C in time dependent CP viola- 
tion for neutral B decay into p'^p~, showing BaBar and Behe results, and the HFAG average. 
Contours are shown at 60.7% confidence level. 

The task at hand is to do a time dependent analysis of the magnitude of the decay amphtudes 
and phases. The analyses of both collaborations allow for p(1450) and p(1700) contributions in 
addition to the /?(770). There are a total of 26 free parameters in the fit. The statistics for 
such an analysis are not overwhelming: BaBar has about 2100 signal events [47], while Belle 
has about 1000 [48]. 

The results for the confidence levels of a found by both collaborations are shown in Figure [T5l 
Note that there is also a mirror solution at a-|-180°. The Belle collaboration uses their 
measured decay rates for p'^ir^ and p^ir^ coupled with isospin relations [49] to 

help constrain a. 

The LHCb experiment expects to be able to significantly improve on the determination of a 
using the pn mode. They expect 14,000 events in for an integrated luminosity of 2 fb~^, with a 
signal to background ratio greater than 1 [50]. It is expected that this amount of data can be 
accumulated in one to two years of normal LHC operation. 

Combining all the data, both the CKM fitter and UT fit groups derive a the confidence level 
plot for a shown in Figure [T6l There is a clear disagreement between the two groups, UT fit 
preferring a solution in the vicinity of 160°, and CKM fitter a value closer to 90°. The CKM 
fitter group believes that this is due to the UT fit group's use of Bayesian statistics which they 
criticize [51]. 
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Figure 13: The probability density, defined as 1 - CL (confidence level) for the angle a as 
determined by measurements oi B ^ pp. 



2.2.5 The Angle 7 



The angle 7 = arg 



. It can be viewed as the phase of the \Vuh\ matrix element, with 
respect to \ Vcb\- Interference measurements are required to determine phases. We can use 
neutral or even charged B decays to measure 7, and there are several ways to do this without 
using any theoretical assumptions, as is the case for a and (3. The first relies on the color 
suppressed tree diagrams shown in Figure [T71 and the second on using the interference in Bg 
mixing. At first glance, the diagrams in Figure [T71 don't appear to have any interfering 
components. However, if we insist that the final state is common to and , then we have 
two diagrams which differ in that one has a b ^ c amplitude and the other a 6 — > n amplitude, 
the relative weak phase is then 7. Explicitly the amplitudes are defined as 



A{B- D'^K-) 



Ab 



(47) 



where rs reflects the amplitude suppression for the b ^ u mode and 6b is the relative phase. 
We have not yet used identical final states for and decays, yet we can see that there will 
be an ambiguity in our result between 8b and 7 that could be resolved by using different 
decay modes. 

There are several suggestions as to different decay modes to use. In the original paper on 
this topic Gronau and Wyler [52] propose using CP eigenstates for the decay, such as 
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Figure 14: Dalitz plot for Monte-Carlo generated p-K ^ -k^tt^-k^ decays. The decays have 

been simulated without any detector effect and the amplitudes for p^vr", p~7r+ and have 
all been generated with equal magnitudes in order to have destructive interferences where the 
p bands overlap. The main overlap regions between the p bands are indicated by the hatched 
areas. Adapted from [47]. 

K'^K~ , vr+vr" etc., combining with charge specific decays and comparing B~ with B'^ decays. 
In the latter, the sign of the strong phase is flipped with respect to the weak phase. In fact 
modes such as D*^ or K*~ can also be used. (When using D*^ there is a difference in 5b of vr 
between jD^ and ir^D^ decay modes [53].) 



It is convenient to define the follow variables: 



A, 



CP± 



Rcp± 



V{B- - 


. r)0(*) ^- 
^ ^CP±^ 


-(*)) _ Y{B+ - 




T{B- - 


^ ^CP±^ 


-(*)) + r(s+ - 





(48) 



These are related to the variables defined in Eq. [57] as 

Rcp± = 1 + ± 2r cos(5b COS7, 
Acp± = ±2r sin 5b sin ^/ Rcp±- 



(49) 



Measurements have been made by the BaBar, Belle and CDF collaborations [54]. These data, 
however, are not statistically powerful enough by themselves to give a useful measurement of 
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Figure 15: The experimental confidence levels for a as determined separately by the Belle 
collaboration (left) and the BaBar collaboration (right). The dotted curve for Belle shows the 
result without using constraints from charged B to pir final states. 
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Figure 16: The experimental confidence levels for a as determined separately by the CKM fitter 
and UT fit groups. 



7. Atwood, Dunietz and Soni suggested using double-Cabibbo suppressed decays as an 
alternative means of generating the interference between decay amplitudes [55]. For example 
the final state B' D^R- can proceed via the tree level diagram in Figure [2]^a) when the 
W~ — > us. Then if we use the doubly-Cabibbo suppressed decay —>■ K^it~ , we get a final 
state that interferes with the D^K~ final state with the — > K'^-k~ , which is Cabibbo 
allowed. The advantage of this method, besides extending the number of useful final states, is 
that both amplitudes are closer to being equal than in the above method, which can generate 
larger interferences. Any doubly-Cabibbo suppressed decay can be used. Similar equations 
exist relating the measured parameters to 7 and the strong phase shift. Measurements have 
been made mostly using the ET^vr^ final state [56]. Again, these attempts do not yet produce 
accurate results. 
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Figure 17: Two diagrams for a charged B decay into a neutral D meson and a charged kaon. 

Thus far, the best method for measuring 7 uses the three-body final state X^vr+Tr^. Since the 
final state is accessible by both and decays interference results. By its very nature a 
three-body state is complicated, and consists of several resonant and non-resonant parts. It is 
typically analyzed by examining the Dalitz plot where two of the possible three sub-masses 
(squared) are plotted versus one another [57]. In this way, the phase space is described as a 
uniform density and thus any resonant structure is clearly visible. For our particular case a 
practical method was suggested by Giri et al. [58], where the decay is analyzed in separate 
regions of the Dalitz plot. Results from this analysis have been reported by BaBar [59] and 
Belle [60]. BaBar finds 7 = (76 ± 22 ± 5 ± 5)°, and Belle finds 7 = (76iJ^ ± 4 ± 5)°. In both 
cases the first error is statistical, the second systematic, and the third refers to uncertainties 
caused by the Dalitz plot model. (There is also a mirror solution at ±180°.) 

It has been shown [61] that measurement of the amplitude magnitudes and phases found in the 
decays of ?/;(3770) (CP it Ta,g){KsT^'^T^~ )D provide useful information that help the narrow 
model error. The CLEO collaboration is working on such an analysis, and preliminary results 
have been reported [62]. (CLEO also is working on incorporating other decay modes.) 

Both the CKM fitter and UT fit groups have formed liklihoods for 7 based on the measured 
results, as shown in Figure [THl In this case CKM fitter and UT fit agree on the general shape 
of the liklihood curve. The CKM fitter plot shows a small disagreement between the Daltiz 
method and a combination of the other two methods. 

2.2.6 The Angle X 

The angle x shown in Figure [T] is the phase that appears in the box diagram for Bs mixing, 
similar to the diagram for B^ mixing shown in Figure [3l but with the d quark replaced by an s 
quark. The analogous mode to B^ J/tpKg in the Bg system is Bg J/tprj. The Feynman 
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Figure 18: The experimental confidence levels for 7 as determined separately by the CKM fitter 
and UT fit groups. 

diagrams are shown in Figure [T9j This is very similar to measuring /3 so x is often called /3<j 




Figure 19: The Feynman diagram for the decay Bg — > J /tpr] or 



Since there are usually two photons present in the r] decay, experiments at hadron colliders, 
which can perform time-dependent studies of Bs mesons, preferentially use the J final 
state. This, unfortunately, introduces another complexity into the problem; as the Bg is 
spinless the total angular momentum of the final state particles must be zero. For a decay into 
a vector and scalar, such as J/ipij, this forces the vector J/tp to be fully longitudinally 
polarized with respect to the decay axis. For a vector-vector final state both angular 
momentum state vectors are either longitudinal (L), both are transverse with linear 
polarization vectors parallel (||) or they are perpendicular (_L) to one another [63]. Another 
way of viewing this is that a spin-0 B decay into two massive vector mesons can form CP even 

■^Note that (j), / -2%, since -2x = argiVtbVtlf / {VcbV'.f , whereas (j), is arg(Mi2/ri2), with arg(Mi2) = 
VtbVtl)^ / {VcbV*^)'^ , and arg(ri2) is a hnear combination of (VcbV's)^ , VcbK'KfcKTs, and (KtKTs)^- 
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states with L=0 or 2, and a CP odd state with L=l. The relative populations in the two CP 
states are determined by strong interactions dynamics, but to study the weak phase here we 
are not particularly interested in the actual amount, unless of course one state dominated. We 
do not expect this to be the case however, since the SU(3) related decay — > J/ipK*^, 
K*^ — > K~^7r~ has a substantial components of both CP states; the PDG quotes gives the 
longitudinal fraction as (80 ± 8 ± 5)% [7]. 

The even and odd CP components can be disentangled by measuring the appropriate angular 
quantities of each event. Following Dighe et al. [64], we can decompose the decay amplitude 
for a Bs as 

A{B, ^ J/V</.) = Aoim^)/E^e*j^^ - A^^efi^l^ - lA^e^ • p/^, (50) 

where ej/^ and are polarization 3- vectors in the J/V' rest frame, p is a unit vector giving 
the direction of the (/> momentum in the J/V' rest frame, and is the energy of the (\) in the 
J/'0 rest frame. We note that the corresponding amplitude for the Bs decay are = A^^ 
A|| = , and A^ = —A±. The amplitudes are normalized so that 

dr{Bs ^ J/ip(l))/dt = \Aof + \A\^f + \A±\^ . (51) 

The (j) meson direction in the J/^|J rest frame defines the x direction. The z direction is 
perpendicular to the decay plane of the K^K~ system , where py{K^) > 0. The decay 
direction of the (.~^ in the J /ip rest frame is described by the angles (^, (p). The angle ijj is that 
formed by the direction with the j;-axis in the (j) rest frame. Figure [20] shows the angles. 



z 




Figure 20: Pictoral description of the decay angles. On the left and <j) defined in the J/ijj rest 
frame and on the right ip defined in the (p rest frame. (From T. Kuhr [65].). 
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The decay width can be written as 



-|- [2 1 Ao p cos^ -0(1 - sin^ 61 cos^ , 

OZTT 



d COS 9 d4> d cos -0 dt 
+ sin^V'{|^|||^(l -sin2 6lsin2 + [^.tl^in^ 61 - Im(A|A_L) sin26'sin0} 

+ ^ sin 2'0{Re(AoA||) sin^ 6 sin 2(?!) + lm.{AlAi_) sin 20 cos cj)}] . 
v 2 

The decay rate for Bg can be found by replacing Aj^ in the above expression with —Aj^. 



(52) 



Another complexity arises from the expectation that the width difference AFs/Fs ~ 15%. This 
complicates the time dependent rate equations. For convenience, setting = (cos9,(p,cosip), 
we have for the decay width for Bg-. 



where 



dtdp 



(X \Ao\^T+h{p) + \Aii\^T+f2{p) 

+ \A^\^T^h{p) + \A^^\\A^\U+U{p) 

+ |Ao||^|||cos(5||)T+/5(p) 

+ \Ao\\A^\V+h{p), 



1 ± cos(2x))e-r^* + (1 T cos(2x))e-r«*] /2, 

ibe~'"* X sin((5_L — 5\\) cos(Amst) 
cos(5_L — (5||) cos(2x) s\n{/S.mst) 
cos((5_L - sin(2x) sinh(AFi/2)] , 



±e~ X [sin((5_L) cos(Amsi) 

— cos((5_l) cos(2x) sin(Amst) 

± cos((5_l) sin(2x) sinh(AFt/2)] . 

/i (p) = 2 cos^ V(l - sin^ ^ cos^ (f)) , 

/2 (p) = sin^ ^ ( 1 — sin^ 9 sin^ <;i!)) , 

fsip) = sm'^i/jsm'^ 9, 

fii^p) = — sin^ V sin(20) sin (f>, 

h{p) = sin(2V')sin2 6'sin(20)/\/2, 

/6(p) = sin(2V')sin(26')cos(/)/A/2. 



(53) 



The quantities 6±_ and are the strong phases of A±_ and A\\ relative to Aq, respectively [66]. 
The expression for Bg mesons can be found by substituting Uj^ — > U- and V+ — >^ V_ . 
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The most interesting quantities to be extracted from the data are x and AF. There are many 
experimental challenges: the angular and lifetime distributions must be corrected for 
experimental acceptances; flavor tagging efficiencies and dilutions must be evaluated; 
backgrounds must be measured. Both the CDF [67] and DO [68] experiments have done this 
complicated analysis. Updated results as of this writing are summarized by the CKM fitter 
derived limits shown in Figure [2T1 




Figure 21: Constraints at 68% confidence level in the {(j)s, AF^) plane. Overlaid are the 
constraints from the CDF and DO measurements, the constraint from AFg = cos (j)sAT^^^ , 
the constraint from the flavor specific Bs lifetime [24], and the overall combination. The SM 
prediction is also given. From [69]. 



The Standard Model allowed region is a very thin vertical band centered near zero at <j)s of 
-0.036ib0.002 (shown in red). The region labeled "all" (green) shows the allowed region at 68% 
confidence level. Although the fit uses several input components besides the CP asymmetry 
measurements in Bs J/tp4', including use of measured total widths introduced via the 
constraint equation AF^ = cos (/)sAFf*^, it is the measurement of ips that dominates the fit 
result. We see that there is a discrepancy that may be as large as 2.7 standard deviations [69]. 
While this is as not yet significant, it is very tantalizing. The LHCb experiment plans to vastly 
improve this measurement [70]. 

It has been pointed out, however, that there is likely an S-wave K~^K~ contribution in the 
region of the (J) that contributes 5-10% of the event rate as estimated using decays [71]. For 
example, analysis of the analogous channel B^ J/tpK*^ reveals about 8% S-wave in the Kit 
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system under the K , and BaBar has used this to extract a value for cos 2/3 thus removing an 
ambiguity in f3 [72]. The S-wave amphtude and phase needs to be added to Eq. \53\ Note that 
the errors wih increase due to the addition of another amphtude and phase. The S-wave can 
manifest itself as a vr+vr", so it is suggested that the decay — > J/'i/'/o(980) be used to 
measure x; here the /o — > ir'^n" [71]; the estimate is that the useful /o rate would be about 
20% of the (j) rate, but the J/ipfo is a CP eigenstate, so an an angular analysis is unnecessary, 
and these events may provide a determination of x with an error comparable to that using 



2.2.7 Measurements of Direct CP Violation in B Ktt Decays 

Time integrated asymmetries of B mesons produced at the T(4S) resonance can only be due to 
direct CP asymmetry, as the mixing generated asymmetry must integrate to zero due to the 
fact that the initial state has J^*^ = 1 . The first evidence for such direct CP violation at 
the greater than four standard deviation level in the K^ir^ final state was given by 
BaBar [73]. The latest BaBar result is [74] 

A.-.. = "'g - ^-"'"> - - = -0.107 ± 0.016^ngS . (54) 

r{B° ^ K~7r+) + r{B0 ^ K+7r~) 

showing a large statistical significance. 

This result was confirmed by the Belle collaboration, but Belle also measured the isospin 
conjugate mode. Consider the two-body decays of B mesons into a kaon and a pion, shown in 
Figure [22l For netural and charged decays, it can proceed via a tree level diagram (a) or a 
Penguin diagram (b). There are two additional decay diagrams allowed for the B~ , the 
color-suppressed tree level diagram (c) and the elusive "Electroweak" Penguin diagram in (d). 
(So named because of the intermediate 7 or Z boson.) Since it is expected that diagrams (c) 
and (d) are small, the direct CP violating asymmetries in both charged and neutral modes 
should be the same. Yet Belle observed [75] 



vr 



0^ 



The Belle data are shown in Figure 

The difference between ^;^-^+ and Axi^nO is not naively expected in the Standard Model and 
Belle suggests that this may be a sign of New Physics. Peskin commented on this 
possibility [76]. 
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2.3 Conclusions from CP Violation Measurements 

All CP violation in the quark sector is proportional to the parameter r] in Eq. [TJ In fact all 
CP asymmetries are proportional to the "Jarlskog Invariant", J = A^X^rj, which represents 
the equal area of all the CKM triangles [77]. Since we know the value of these three numbers, 
we do know the amount of CP violation we can expect, even without making the 
measurements. We also can estimate the amount oi CP violation necessary using cosmology. 
To reproduce the observed baryon to entropy ratio requires many orders of magnitude more 
CP violation than thus far found in heavy quark decays [78]. Thus we believe there are new 
sources of CP violation that have not yet been found. 



3 The CKM Parameter \Vcb\ 

There are two experimental methods to determine \ Vcb\'- the exclusive method, where \ Vcb\ is 
extracted by studying the exclusive B D^*^£~i? decay process; and the inclusive method, 
which uses the semileptonic decay width of 6-hadron decays. In both methods, the extraction 
of \Vcb\ is systematics limited and the dominant errors are from theory. The inclusive and 
exclusive determinations of \ Vcb\ rely on different theoretical calculations of the hadronic 
matrix element needed to extract it from measured quantities, and make use of different 
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Figure 23: The beam constrained mass distributions from Belle in units of GeV for the four 
different indicated final states. The dotted curves indicate peaking backgrounds, the dot-dashed 
curves the signal and the solid curves the sum. The differences in event numbers between the 
charge conjugate modes are apparent. 



techniques which, to a large extent, have uncorrelated experimental uncertainties. Thus, the 
comparison between inclusive and exclusive decays allows us to test our understanding of 
hadronic effects in semileptonic decays. The latest determinations differ by more than 2cr, with 
the inclusive method having a stated error half of the size of the exclusive one. 



3.0.1 Beauty Quark Mass Definitions 

Due to confinement and the non-perturbative aspect of the strong interaction, the concept of 
the quark masses cannot be tied to an intuitive picture of the rest mass of a particle, as for 
leptons. Rather, quark masses must be considered as couplings of the SM Lagrangian that 
have to be determined from processes that depend on them. As such the 6-quark mass (mf,) is 
a scheme-dependent, renormalised quantity. 

In principle, any renormalisation scheme or definition of quark masses is possible. In the 
framework of QCD perturbation theory the difference between two mass schemes can be 
determined as a series in powers of a^. Therefore, higher-order terms in the perturbative 
expansion of a quantity that depends on quark masses are affected by the particular scheme 



33 



employed. There are schemes that are more appropriate and more convenient for some 
purposes than others. Here we examine the main quark mass definitions commonly used in the 
description of B decays. 



• Pole mass: The pole mass definition is gauge-invariant and infrared-safe [79] to all 
orders in perturbation theory and has been used as the standard mass definition of many 
perturbativc computations in the past. By construction, it is directly related to the 
concept of the mass of a free quark. The presence of a renormalon ambiguity [80] makes 
the numerical value of the pole mass an ordcr-dcpcndcnt quantity, leading to large 
perturbative corrections for Heavy Quark Effective Theory (HQET) parameters (see 
below for a discussion of HQET). These shortcomings are avoided by using quark mass 
definitions that reduce the infrared sensitivity by removing the Aqcd renormalon of the 
pole mass. Such quark mass definitions are generically called "short-distance" masses. 

• MS mass: The most common short-distance mass definition is the MS mass rnb(/Li) [81], 
where the scale fi is typically chosen to be the order of the characteristic energy scale of 
the process. In the MS scheme the subtracted divergencies do not contain any infrared 
sensitive terms and the MS mass is only sensitive to scales larger than mi,. The MS mass 
arises naturally in processes where the 6-quark is far off-shell, but it is less adequate 
when the 6-quark has non-relativistic energies. 

• Kinetic mass: The shortcomings of the pole and the MS mass in describing 
non-relativistic 6-quarks can be resolved by so-called threshold masses [82] , that are free 
of an ambiguity of order Aqcd and are defined through subtractions that contain 
universal contributions for the dynamics of non-relativistic quarks. Since the 
subtractions are not unique, an arbitrary number of threshold masses can be 
constructed. The kinetic mass is defined as [83,84]: 



where /^kin is the nominal kinetic mass renormalisation scale. For /^kin the kinetic 
mass reduces to the pole mass. 

• IS' mass: The kinetic mass depends on an explicit subtraction scale to remove the 
universal infrared sensitive contributions associated with the non-relativistic 6-quark 
dynamics. The IS mass [85] achieves the same task without a factorisation scale, since it 
is directly related to a physical quantity. The 6-quark IS mass is defined as half of the 
perturbative contribution to the mass of the T{Si) in the limit 
mj, » rribV » mbv"^ » Aqcd- 



"i6,km(Mkiii) — "Ifo.pole ~ [A (//kin)] pert 




+ ..., 



(56) 



2m6,kin(//kin) 



- pert 
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Table 2: List of rrn, determinations converted into the MS mass scheme. 



mh{mi,) GeV 



Method 



4.243 ± 0.042 
4.25 ±0.02 ±0.11 
4.346 ± 0.07 
4.164 ±0.025 
4.17 ±0.05 



From B Xs^f and B XJu fit [86] 



Lattice UKQCD [87] 
T(IS') NNNLO [88] 



low-moment sum rules NNNLO [89] 
T sum rules NNLO [90] 



A list of &-quark mass determinations, converted into the MS mass scheme is shown in Tabled 
3.1 Determination Based on Exclusive Semileptonic B Decays 

The exclusive |V^f,| determination is obtained by studying the decays B D*(.v and B Diu, 
where i denotes either an electron or a muon. The exclusive measurements of a single hadronic 
final state, e.g. the ground state D or D* , restrict the dynamics of the process. The remaining 
degrees of freedom, usually connected to different helicity states of the charmed hadron, can be 
expressed in terms of form factors, depending on the invariant mass of the lepton-j/ pair, q'^. 
The shapes of those form factors are unknown but can be measured. However, the overall 
normalization of these functions needs to be determined from theoretical calculations. 

Isgur and Wise formulated a theoretical breakthrough in the late 1980's. They found that in 
the limit of an infinitely heavy quark masses QCD possess additional flavor and spin 
symmetries. They showed that since the heavier b and c quarks have masses much heavier 
than the scale of the QCD coupling constant, that they are heavy enough to posses this 
symmetry but that corrections for the fact that the quark mass was not infinite had to be 
made. They showed there was a systematic method of making these corrections by expanding 
a series in terms of the inverse quark mass [91]. This theory is known as Heavy Quark 
EflFective Theory (HQET). 

When studying b quark decays into c quarks, it is convenient to view the process in 
four-velocity transfer (w) space as opposed to four-momentum transfer space, because at 
maximum four-velocity transfer, where w equals one, the form-factor in lowest order is unity, 
i.e. the b transforms directly into a c quark without any velocity change. The value of \Vcb\ can 
be extracted by studying the decay rate for the process B^ D*^i^V£ as a function of the 
recoil kinematics of the D*'^ meson. Specifically, HQET predicts that: 



dT{B ^ DHv) _ Gl\Vcb? 



(57) 
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where w is the product of the four- velocities of the D*+ and the mesons, related to q^. 
)C{w) is a known phase space factor and the form factor J^{w) is generally expressed as the 
product of a normalization factor and a shape function described by three form factors, 
constrained by dispersion relations [92]: 



hAi (w) 



w - rf (1 - R2{w)f + 2(1 - 2wr - r^) ( 1 - - — \r^{w) 



(l-r)2 

1 - Sp^z + (53/9^ - lh)z^ - (231/9^ - 91)z^, 



w + l 



,(58) 



/iAi(l) 

R^{w) = i?i(l) - 0.12(u; - 1) +0.05(u' - 1)^ 
R^{w) = i?2(l) - 0.11(u; - 1) - 0.06(u' - 1)2, 

where r = mo* /i^b, z = ^^E^^; The linear slope of the form- factor is given by the 
parameter p^, and must be determined from the data. In the infinite mass limit, 
J-'{w = 1) = hAiiu! = 1) = 1; for finite quark masses, non-pcrturbativc effects can be expressed 
in powers of l/rriQ. There are several different corrections to the infinite mass value J-{1) = 1: 

^(1) = r/QEDT/A [l + h/ml + ••■] (59) 

Note that the first term in the non-perturbative expansion in powers of l/mg vanishes [93]. 

QED corrections up to leading logarithmic order give rjA = 0.960 it 0.007. Different estimates 

of the l/rriQ corrections, involving terms proportional to 1/m^, l/m^ and l/{mbmc), have 

been performed in a quark model with QCD sum rules, and, more recently, with an HQET 

based lattice gauge calculation. The best estimate comes from lattice QCD, 

hAiil) = J^(l) = 0.921 ± 0.013 ± 0.020 [94]. This result does not include a 0.7% QED 

correction. 

Since the phase-space factor K.{w) tends to zero as w ^ 1, the decay rate vanishes and the 
accuracy of the \Vcb\ value extracted with this method depends upon experimental and 
theoretical uncertainties in the extrapolation. Experiments determine the product |J^(1) • K;6p 
by fitting the measured dT/dw distribution. 

This decay has been analyzed by CLEO, BaBar and Belle using B mesons from the T(4S) 
decay, and by ALEPH, DELPHI, and OPAL at the Z° center of mass energy. Experiments 
that exploit the T(4S) have the advantage that w resolution is quite good. However, they 
suffer from lower statistics near w = 1 in the decay B D*^l!.u due to the lower 
reconstruction cflicicncy of the slow vr^. On the other hand, the decay B D*^lu is not 
affected by this problem [95]. In addition, kinematic constraints enable these experiments to 
identify the final state including D* without large contamination from the poorly known 
portion of semileptonic decays with a hadronic system recoiling against the lepton-i^ pair with 
masses higher than the D and £)*, commonly identified as '£)**'. B -factories and CLEO fit for 
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Table 3: Experimental results based on B ^ D*iv after the correction to common inputs and 
the world average as obtained by HFAG [99]. All numbers are corrected to use Ri and R2 from 
a global fit to CLEO and BaBar data. The newest BELLE result is not included in the world 
average. 



experiment 


H'^)\v,k\ (xio-^) 




Corrstat 


Branching fraction(%) 


ALEPH [100] 


31.6±1.8±1.3 


0.50±0.20±0.09 


94% 


5.44±0.25±0.20 


BaBar [101] 


33.9±0.3±1.1 


1.18±0.05±0.03 


27% 


4.53±0.04±0.13 


BaBar [102] 


34.9±0.8±1.4 


1.11±0.06±0.08 


90% 


5.40±0.16±0.25 


BaBar [103] 


35.7±0.2± 1.2 


1.20±0.02±0.07 


38% 




BELLE [97] 


34.7± 0.2± 1.0 


1.16± 0.04± 0.03 


91% 


t 


BELLE [104] 


34.7± 0.2±1.0 


1.16±0.04±0.03 


91% 


4.75±0.25± 0.19 


CLEO [105] 


41.3±1.3±1.8 


1.37±0.08±0.18 


91% 


6.02± 0.19±0.20 


DELPHI(excl) [100] 


36.3±1.8±1.9 


1.04±0.14±0.15 


89% 


5.53±0.19±0.34 


DELPHI(part rec) [100] 


35.8±1.4±2.3 


1.18±0.13± 0.25 


94% 


5.00±0.15±0.18 


OPAL(excl) [100] 


36.9±1.6±1.5 


1.24±0.20±0.14 


77% 


5.17± 0.20±0.38 


OPAL(par rec) [100] 


37.6±1.2±2.4 


1.14±0.13±0.27 


94% 


5.63±0.27± 0.43 


world average 


35.4 ±0.5 


1.16 ±0.05 


20% 


5.05 ± 0.10 



f - Not used in the World Average 



the signal and background components in the distribution of the cosine of the angle between 
the direction of the B and the direction of the D*l system. At LEP, B mesons are produced 
with a large variable momentum (about 30 GeV on average), giving a relatively poor w 
resolution and limited physics background rejection capabilities. By contrast, LEP 
experiments benefit from an efficiency that is only mildly dependent upon w. 

LEP experiments extracted \V(^\ by performing a two-parameter fit, for .F(l)|V^b| ^-^d the slope 

. The first measurements of both ratio Ri and i?2, and , were made by the CLEO 
collaboration [96]. Belle [97] and BaBar [98] improved upon these measurements using the 
B^ D*^e~h'(. decay. They determined R2, and using an unbinned maximum 
likelihood fit to the full decay distribution. BaBar and CLEO results are combined to give: 
Ri = 1.396 ± 0.060 ± 0.035 ± 0.027, R2 = 0.885 ± 0.040 ± 0.022 ± 0.013, and 

= 1.145 ± 0.059 ± 0.030 ± 0.035. The stated uncertainties are the statistical from the data 
and systematic uncertainty, respectively. 

Table [3] summarizes the available data. Values of .F(1)|V^;,| from different experiments can be 
combined if they are extracted using the same J'{w) parametrization. All measurements 
included in the exclusive \ Vch\ world average relying on the above form factor ratios i?i and i?2- 
Using the J^(l)|V;fe| world average in Tableland /iai(1) = 0.921 ± 0.013 ± 0.020 [94], we find 
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Table 4: Experimental results from B — > Div using corrections to the common inputs and world 
averages. 



experiment 


g(l)|Veb|(xlO-^) 


Pi 


Corrstat 


Branching fraction(%) 


ALEPH [100] 


38.1±11.8±6.1 


0.91±0.98± 0.36 


98% 


2.25 ±0.18 ±0.36 


BaBar [107] 


42.3±1.9±1.0 


1.20±0.09±0.04 


95% 




BaBar [108] 


43.8±0.8±2.1 


1.22±0.04±0.06 


63% 




BaBar [109] 








2.20 ±0.11 ±0.12 


BELLE [110] 


40.7±4.4±5.1 


1.12±0.22±0.14 


96% 


2.09 ±0.12 ±0.39 


CLEO [111] 


44.5±5.9±3.4 


1.27±0.25±0.14 


95% 


2.10 ±0.13 ±0.15 


World average 


42.4 ± 0.7 ± 1.4 


1.19 ±0.04 ±0.04 


96% 


2.16 ±0.12 



\Vcb\ = (38.2 ± O.bexp ± l-Otheo) X 10-3, 
where the dominant error is theoretical, and it will be difficult to improve upon. 

The study of the decay B Dlv poses new challenges from the experimental point of view. 
The differential decay rate for B Dlv can be expressed as [106] 

- = ^l^/c. HaH^ (60) 

where ICx>{w) is the phase space factor and the form factor G{w) is generally expressed as the 
product of a normalization factor ^(1) and a function, g£){w), constrained by dispersion 
relations [92]. The strategy to extract ^(l)|Vcb| is identical to that used for the B — > D*ii' 
decay. However, dTu/dw is more heavily suppressed near w = 1 than dL^)* /dw due to the 
helicity mismatch between initial and final states. Moreover, this channel is much more 
challenging to isolate from the dominant background B D^iv as well as from fake D-l 
combinations. Table [4] shows the results of two-dimensional fits to |T4fe|^(l) and for different 
experiments and the world average. 

In the limit of infinite quark masses, G{w = 1) coincides with the Isgur-Wise function [112]. In 
this case there is no suppression of l/mg corrections and QCD effects on ^(1) are calculated 
with less accuracy than .F(l). Corrections to this prediction have recently been calculated with 
improved precision, based on unquenched lattice QCD [113], specifically 
g{l) = 1.074 ± 0.018 ± 0.016. Using this result we get 

\Vcb\ = (39.5 ± 1.4e:,p ± 0.9theo) X 

consistent with the value extracted from B D*^^ decay, but with an experimental 
uncertainty about twice as large. 
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BaBar recently has also studied the differential decay widths for the decays B D lu and 
B- D*°£u to extract the ratio ^(1)/J^(1) = 1.23 ± 0.09 [114], compatible with the lattice 
theory prediction of 1.16 it 0.04 . 

3.2 B D**lu Decays 

It is important to understand the composition of the inclusive B semileptonic decay rate in 
terms of exclusive final states for use in semileptonic B decay analyses. The B D^*^ii' 
decays are well measured, but a sizeable fraction of semileptonic B decay are to D**lv. The 
D** resonances have larger masses than the D* and not well studied^] The measurements of 
these branching fractions require a good understanding of the different D^*^mr systems, which 
can be either resonant or non-resonant, each with characteristic decay properties. 

There are four orbitally excited states with L=l. They can be grouped in two pairs according 
to the value of the spin on the light system, j =Libl/2 (L=l). States with j = 3/2 can have 

= l"''and 2"*" state. The 1"*" state decays only through D*7r, and the 2"^ through Dtt or D*Tr. 
Parity and angular momentum conservation imply that in the 2+ the D* and vr are in a D 
wave but allow both S and D waves in the 1"*" state. However, if the heavy quark spin is 
assumed to decouple, conservation of j = 3/2 forbids S waves even in the l"*" state. A large 
D-wave component and the fact that the masses of these states are not far from threshold 
imply that the j = 3/2 states are narrow. These states have been observed with a typical 
width of 20 MeV/c^. On the contrary, j = 1/2 states can have = O^and 1+, so they are 
expected to decay mainly through an S wave and manifest as broad resonances, with typical 
widths of several hundred MeV/c^. 

The ALEPH [100], CLEO [115], DELPHI [100], and DO [116] experiments have reported 
evidence of the narrow resonant states {Di and D2) in semileptonic decays, whereas more 
recent measurements by the BaBar [117] and Belle [118] experiments provide semi-inclusive 
measurements to D^*^ttIu final states [119]. 

The differences between the measured inclusive semileptonic branching fraction and the sum of 
all exclusive B semileptonic measurements for Belle, BaBar and World averages are given in 
Table [5] for B^ and B~ decays. In the case where multiple measurements exist, only the most 
precise measurements have been used in the BaBar and Belle columns, i.e. no attempt at an 
average is made. In all cases the sum of the exclusive components does not saturate the B 
semileptonic rate. 

^ D** refers to the resonant states, Dq, Dl, D2 and Di. 
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Table 5: Inclusive versus sum of exclusive measured B semileptonic branching fractions (%). 



B{%) 


BaBar 
B° B- 


Belle 

B° B- 


World Average 
B° B- 


B Div 
B ^ D*lu 
B Dniv 
B D'ttIv 


2.20 ± 0.16 2.30 ± 0.10 
4.53 ± 0.14 5.37 ± 0.21 
0.42 ± 0.09 0.42 ± 0.07 
0.48 ± 0.09 0.59 ± 0.06 


2.09 ± 0.16 
4.42 ± 0.25 

0.43 ± 0.09 0.42 ± 0.06 
0.57 ± 0.22 0.68 ± 0.11 


2.16 ± 0.12 2.32 ± 0.09 
5.05 ± 0.10 5.66 ± 0.18 
0.43 ± 0.06 0.42 ± 0.05 
0.49 ± 0.08 0.61 ± 0.05 


S(Exc.) 
Inc. 


7.63 ± 0.25 8.68 ± 0.25 
10.14± 0.43 10.90 ± 0.47 


7.51 ± 0.73 

10.46 ± 0.38 11.17 ± 0.38 


8.13 ± 0.19 9.01 ± 0.21 
10.33 ± 0.28 10.99 ± 0.28 


Inc. - S(Exc.) 


2.51 ± 0.50 2.22 ± 0.53 


2.95 ± 0.82 


2.20 ± 0.34 1.98 ± 0.35 



All measured rates for the D** narrow states are in good agreement. Experimental results 
seem to point towards a larger rate for broader states. If it is due mainly to D'l and Dq decay 
channels, these results disagree with the prediction of QCD sum rules. However, Belle set an 
upper limit for the D[ channel below the rate measured by BaBar and DELPHI. More 
measurements need to be performed to elucidate this puzzle. 

3.3 Determination Based on Inclusive Semileptonic B Decays 

Inclusive determinations of \Vcb\ are obtained using combined fits to inclusive B decay 
distributions [120,121]. These determinations are based on calculations of the semileptonic 
decay rate in the frameworks of the Operator Product Expansion (OPE) [122] and 
HQET [120, 123]. They predict the semileptonic decay rate in terms of jV^bj, the 6-quark mass 
rrih, and non-perturbative matrix elements. The spectator model decay rate is the leading term 
in a well-defined expansion controlled by the parameter Aqcd/''^6 [123-126] with 
non-perturbative corrections arising to order l/rn^. The key issue in this approach is the 
ability to separate perturbative and non-perturbative corrections (expressed in powers of as). 
Thus the accuracy of an inclusive determination of \Vcb\ is limited by our knowledge of the 
heavy quark nonperturbative matrix elements and the b quark mass. It also relies upon the 
assumption of local quark hadron duality. This is the statement that a hadronic matrix 
element can be related pointwise to a theoretical expression formulated in terms of quark and 
gluon variables. 

Perturbative and non-perturbative corrections depend on the mf, definition, i.e. the 
non-perturbative expansion scheme, as well as the non-perturbative matrix elements that enter 
the expansion. In order to determine these parameters. Heavy Quark Expansions 
(HQE) [120,126,127] express the semileptonic decay width Psl, moments of the lepton energy 
and hadron mass spectra in B ^ Xcif decays in terms of the running kinetic quark masses 
m^™ and m^'" as well as the 6-quark mass in the IS expansion scheme. These schemes 
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should ultimately yield consistent results for \Vcb\- The precision of the 6-quark mass is also 
important for \Vub\- 

The shape of the lepton spectrum and of the hadronic mass spectrum provide constraints on 
the heavy quark expansion, which allows for the calculation of the properties of -B — > Xcli^ 
transitions. So far, measurements of the hadronic mass distribution and the leptonic spectrum 
have been made by BaBar [128], Belle [129], CLEO [130,131], DELPHI [132]. CDF [133] 
provides only the measurement of the hadronic mass spectrum with a lepton momentum cut of 
0.6 GeV in the B rest frame. 



The inclusive semileptonic width can be expressed as 

where p = m^/ml, and /u^, Hq, pd and p^s are non-perturbative matrix elements of local 
operators in HQET. To make use of equations such as Eq. [611 the values of the 
non-perturbative expansion parameters must be determined experimentally. Although some of 
these parameters are related to the mass splitting of pseudoscalar and vector meson states, 
most non-perturbative parameters are not so easily obtained. Measurements of the moments of 
different kinematic distributions in inclusive S-decays are used to gain access to these 
parameters. The first moment of a distribution is given by: 

(Ml) = j Mi{x)dx—, (62) 

corresponding to the mean. Subsequent (central) moments are calculated around the first 
moment, 

(M„) = ((f-Mi)"), (63) 
corresponding to a distribution's width, kurtosis and so on. 

A B ^ Xclv decay observable calculated with the OPE is a double expansion in terms of the 
strong coupling as{mb) and the ratio Aqcd/^ti;,. Observables are typically calculated with the 
cuts used in the experimental determination for background suppression, and enhanced 
sensitivity such as the lepton energy cut. The spectral moments are defined as: 

<EfKiM\)'> = ^ 4,,, m j ^.E. j iM, ^^J^^^^^ E'iEf(M\)' (64) 
= /o[n,E,] + /i[n,E,]A 
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where 



, _ Glml\V,,\' 
° 1927r3 



(65) 



The measurement of \Vcb\ from inclusive decays requires that these decays be adequately 
described by the OPE formalism. The motivation of the moment approach is to exploit the 
degree of experimental and theoretical understanding of each moment and for different Ef^^^, 
directly examining the dependence of the various coefficient functions /j on these terms. The 
possibility of deviations from the OPE predictions due to quark hadron duality violations have 
been raised [134]. To compare the OPE predictions with data, one also has to define how 
uncertainties from corrections are estimated. These uncertainties are hard to quantify 

reliably, as the only information on the size of the matrix elements of the dimension six 
operators comes from dimensional analysis. The only method to check the reliability of \Vcb\ 
extraction from inclusive semileptonic decays is how well the OPE fits the data. 

To compare with theoretical predictions, the moments are measured with a well defined cut on 
the lepton momentum in the B rest frame. The measured hadronic mass distribution and 
lepton energy spectrum are affected by a variety of experimental factors such as detector 
resolution, accessible phase space, radiative effects. It is particularly important to measure the 
largest fraction of the accessible phase space in order to reduce both theoretical and 
experimental uncertainties. Each experiment has focused on lowering the lepton energy cut. 

The hadronic mass spectrum in B ^ Xclu decays can be split into three contributions 
corresponding to D, D*, and D**, where D** here stands for any neutral charmed state, 
resonant or not, other than D and D*. Belle [129], BaBar [128] and CLEO [130] explored the 

moments of the hadronic mass spectrum as a function of the lepton momentum cuts. 
CLEO performs a fit for the contributions of signal and background to the full 
three-dimensional differential decay rate distribution as a function of the reconstructed 
quantities q^, M^, cos^^. Belle [135] and BaBar [128] use a sample where one of the B mesons, 
produced in pairs from T(45) decays, is fully reconstructed and the signal side is tagged by a 
well identified lepton. The 4-momentum px of the hadronic system X, recoiling against the 
lepton and neutrino, is determined by summing the 4-momenta of the remaining charged tracks 
and unmatched clusters. Belle reconstructs the full hadronic mass spectrum and measures the 
first, second central and second non-central moments of the unfolded Mj- spectrum in 
B -^■ Xcii', for lepton energy thresholds, -Emim varying from 0.7 to 1.9 GeV [129]. BaBar 
extracts the moments from the measured distributions using a calibration curve derived from 
Monte Carlo data, with a minimum momentum for the electron in the B meson rest frame of 
0.9 GcV. The latest BaBar analysis [129] measures the first, second central and second 
non-central moments of the spectrum for -Bmin from 0.9 to 1.9 GeV. The main systematic 
errors originate from background estimation, unfolding and signal model dependence. 
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DELPHI follows a different approach in extracting the moments, measuring the invariant mass 
distribution of the D** component only and fixing the D and D* components. DELPHI 
measures the first moment with respect to the spin averaged mass of D and D*. At LEP 
6-quarks were created with an energy of approximately 30 GeV allowing the measurement of 
the hadronic mass moments without a cut on the lepton energy [132]. 

The shape of the lepton spectrum provides further constraints on the OPE. These 
measurements are sensitive to higher order OPE parameters and are considerably more precise 
experimentally. Moments of the lepton momentum with a cut P£ > 1.0 GeV/c have been 
measured by the CLEO collaboration [136]. BaBar [128] extract up to the third moment of 
this distribution, using a low momentum cut of pi > 0.6 GeV/c. Both BaBar and CLEO use 
dilepton samples. The most recent measurement of the electron energy spectrum is from 
Belle [135]. Events are selected by fully reconstructing one of the B mesons, produced in pairs 
from T(4S') decays and it determines the true electron energy spectrum by unfolding [137] the 
measured spectrum in the B meson rest frame. Belle measures B^ and B^ weighted average 
partial branching fractions B{B ^ Xcll.y)Ei>E^i^ and the first four moments of the electron 
energy spectrum m. B ^ Xceu, for i^min from 0.4 to 2.0 GeV [138]. All the lepton moment 
measurements are consistent with theory and with the moment of the hadronic and 6 — > 57 
photon energy spectrum. Hadronic and lepton energy measurements are consistent within 
their errors. When compared with theory there is no sign of inconsistencies. The B-factories 
provide the most precise HQE parameter estimates. 



3.3.1 HQE Parameters 

Using the moment measurements described above, it is possible to determine the CKM matrix 
element \ Vcb\ and HQE parameters by performing global fits in the kinetic and IS 6-quark mass 
schemes [139]. The photon energy moments m. B ^ decays [140] are also included in 
order to constrain the 6-quark mass more precisely. Measurements that are not matched by 
theoretical predictions and those with high cutoff energies are excluded (i.e. semileptonic 
moments with -Emin > 1-5 GeV and photon energy moments with -Emin > 2 GeV). The results 
are preliminary. 

The inclusive spectral moments of ^ Xciv decays have been derived in the IS scheme up to 
0{l/ml) [120]. The theoretical expressions for the truncated moments are given in terms of 
HQE parameters with coefficients determined by theory, as functions of i^min- The 
non-perturbative corrections are parametrized in terms of the following non-perturbative 
parameters: A at 0{mb), Xi and A2 at 0(l/m^), and ti, T2, T3, T4, pi and p2 at 0{l/m\). In 
Table [6] one finds the following results for the fit parameters [139]. The first error is from the 
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Table 6: Experimental results using corrections to the common inputs and world averages for 
the IS scheme. 





l^cfcKxio-^) 


mf (GeV) 


Ai(GeV^) 


X^/n.d.f. 




~X^7 41.81 ±0.34fit±0.08rB 


4.700 ±0.030 


-0.315 ±0.026 


24.7/56 




42.03 ± 0.42fit ± O.OSrs 


4.656 ±0.060 


-0.343 ±0.046 


19.0/45 




Figure 24: Fit results for \Vch\ and to i? ^ Xcif data only (larger blue ellipse) and 
B Xclv and B Xs^ data combined (smaller red ellipse) for the IS" scheme. The ellipses 
are for Ax^ = 1 [141]. 



fit including experimental and theory errors, and the second error (on \Vcb\ only) is due to the 
uncertainty on the average B lifetime. If the same fit is performed to all measured moments of 
inclusive distributions va. B ^ X^lf and B ^ decays, the \ Vcb\ and values obtained are 
in Table [6l The fit results for \ Vch\ and to B ^ Xcly data only and B Xciv and 
B Xgj data combined are displayed in Fig. [241 

Spectral moments of B ^ Xcd-v decays have been derived up to 0{l/m\) in the kinetic 
scheme [126]. The theoretical expressions used in the fit contain improved calculations of the 
perturbative corrections to the lepton energy moments [84] and account for the 
Erain dependence of the perturbative corrections to the hadronic mass moments [142]. For the 
B moments, the (biased) OPE prediction and the bias correction have been 

calculated [127]. All these expressions depend on the h- and c-quark masses mi,{^) and mdp), 
the non-perturbative parameters, defined at the scale ^ = 1 GeV: )U^(/i) and fiQifi) (©(l/m^)), 
p\){p) and Pisif^) (C'(l/''7^5)), and a^. The CKM element \Vcb\ is a free parameter in the fit. 
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Table 7: Experimental results using corrections to the common inputs and world averages for 
the kinetic scheme. The first error is from the fit (experimental error, non-perturbative and bias 
corrections), and the second error is a 1.5% uncertainty from the theoretical expression for the 
semileptonic width [123]. 

\Vcb\{xW-'') rrib (GeV) l^ljGeY'') xVn.d.f. 

Xciiy + Xsj 41.67 ± 0.47/it ibO.OSrs ±0.58t/, 4.601 ± 0.034 0.440 ± 0.040 29.7/57 
XJu + Xsj] 40.85 ±0.68/it± 0.08^5 ±0.57t/, 4.605 ± 0.031 0.312 ± 0.060 54.2/46 
XJu 41A8 ±0A7 fit ± O.OSrB ± 0.58th 4.659 ifc 0.049 0.428 ifc 0.044 24.1/46 

f - New theoretical correlations are applied [143]. 

related to the semileptonic width T{B Xciv) [123]. The fit results for various inputs and 
correlations are shown in Table [7] . 

All the measured moments of inclusive distributions in — > Xciv and B ^ sj decays are used 
in a fit to extract \Vcb\ and the b and c quark masses. The \ Vcb\ and values obtained are 
listed in the Tabled The default fit also gives rUc = 1.16 ± 0.05 GeV, and /i^ = 0.27 ± 0.04 
GeV^. The errors are experimental and theoretical (HQE and Fsl) respectively. In this fit the 
following variations were considered when in the extraction of the HQ parameters, ±20 MeV 
for the b and c quark masses, ±20% for 1^1 and /i^, ±30% for the 3'^'=' order non perturbative 
terms and 02 = 0.22 ± 0.04 for the perturbative corrections. The bias corrections uncertainties 
for i? — > 57 were varied by the full amount of their value. 

There are open issues relating to the global fits. First of all, the x^/ n.d.f. are very small, 
pointing to an underestimate in the theoretical correlations. In a recent study [143], the 
theoretical correlations used in the fit were scrutinized, and new correlation coefficients were 
derived from the theory expressions using a "toy Monte Carlo" approach, showing that the 
theoretical correlations were largely underestimated. The result of this new fit is shown in 
Table [71 The second issue is related the size of the theoretical error. Recently the NNLO full 
two-loop calculations become available [144] . In the Kinetic scheme NNLO calculations include 
an estimate of the non-BLM terms and lead to a roughly 0.6 % reduction of the \Vcb\ value 
—0.25 X lO^'^. In the IS scheme the shift on \Vcb\ is of about —0.14 x 10"'^ [145]. From the new 
power corrections at NLO we expect the chromo-magnetic corrections to be more important as 
the tree level corrections are more important, and a change of about 20-30% in the extracted 
value of //^ in the pole expansion (may be less with the other schemes) [146]. HQE has been 
carried out up to 1/m^ and the effects are expected to be of the order J^^^F/F ^ 0.25% [147]. 
All these newly calculated results can be used to scrutinize the earlier error estimates. In the 
kinetic scheme, the full NNLO value for Ap^rt = 0.919, which is in good agreement with 
Apert = F(S XJv)/T{B XJu)tree = 0.908 ± 0.009. In the IS scheme the estimated 
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uncertainty from the non-BLM two-loop is half of the BLM part, equivalent to 1.5% of the 
tree-level rate, more than 3 times the actual correction. The new results are not implemented 
in the global fits as they are available mainly in the pole-mass scheme. 



3.4 Outlook 

The error on the inclusive and exclusive determination of Vcb is limited by theory. The two 
method have a 2a disagreement, even when results in the same experiment are compared. 

For exclusive determinations, the experimental determinations m B ^ D*ll.u from different 
experiments are not in agreement. Improving the statistical error for the determination from 
B — > Div will help in elucidating the origin of this discrepancy. 

In the inclusive method, errors of less than 2% are quoted. However, the latest theoretical 
results and the introduction of better correlation between theoretical error show shifts in the 
central value larger than the quoted fit error. The situation needs to be reevaluated when all 
the new calculations and corrections are implemented in the fit. Another puzzling result of the 
global fit to i? — > Xclv moments only is the value of m^'^, which seems to be in disagreement 
with other determinations, such as the ones in Table Wn 



4 The CKM parameter |Kfe| 

The parameter \Vub\ determines one of the sides of the unitarity triangle, and thus affects one 
of the crucial tests of the Yukawa sector of the Standard Model. Also in this case, there are 
two general methods to determine this parameter, using B meson semileptonic decays. The 
first approach relies on the determination of branching fractions and form factor 
determinations of exclusive semileptonic decays, such as ^ niUi. The relationship between 
experimental measurements and \Vub\ requires a theoretical prediction of the hadronic form 
factors governing these decays. The complementary approach relies on measurements of 
inclusive properties of B meson semileptonic decays. In this case, the hadronic matrix element 
is evaluated in the context of the OPE. 

Both methods pose challenges to both experimenters and theorists. The branching fractions 
are small, and a substantial background induced by the dominant h ciui needs to be 
suppressed or accounted for. The large data samples accumulated at the two b-factories, Belle 

''For a detailed comparison between m^'" and different mass determinations see A. ffoang, talk at Joint 
Workshop on \Vub\ and |Vc6| at the B- Factories Heidelberg, December f4-f6, 2007. 
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and BaBar, have made possible the development of new experimental techniques that have 
reduced the experimental errors substantially. In parallel, theorists devoted considerable efforts 
to devise measurable quantities for which reliable predictions could be produced, and to 
improve the accuracy of the predictions through refinements in the calculation. Although the 
precision of the stated errors improved, so did the difference in central values between the 
inclusive and exclusive estimates of \Vuh\, at least in most determinations. Possible 
interpretations of this discrepancy will be discussed. 

4.1 Determinations Based on Exclusive B Semileptonic Decays 

The decay B ttIui is the simplest to interpret, as it is affected by a single form factor. The 
differential decay width is given by 

where Gp is the Fermi constant, A(g^) = (g^ + m\ — m^)^ — Am^ml., and f+{(f) is the 
relevant form factor. 

The first experimental challenge is to reconstruct this exclusive channel without significant 
background from the dominant charm semileptonic decays, and the additional background 
component from other h u£i> transitions. The advantage of e+e" B-factories is that the B 
decaying semileptonically originates from the decay e'^e~ BB. Thus, if the companion B is 
fully reconstructed, the 4-momentum can be inferred from the rest of the event. 

CLEO pioneered this approach by reconstructing the i' from the missing energy 

{Ejniss = 2£^s — TiiEj), and momentum {pmiss = ^iPi) in the event; in these definitions the 

index i runs over the well reconstructed tracks and photons, and cuts are applied to enhance 

the probability that the only missing particle in the event is the i>. This approach allows the 

application of a low momentum cut of 1.5 GeV on the lepton for B irive and 2 GeV for 

B p£i>e. Averaging the results from these two exclusive channels, they obtain 

\Vub\ = (3.3 ± 0.21^:1 ± 0.7) X 10-3. Using their relatively modest fuh data set (16 fb^^) at the 

center-of-mass energy of the T(45') and considering only the B vr^P^ channel, they get the 

branching fraction B{B Tr+lue) = (1.37 ± 0.15 it 0.11) x 10"^, and 

\Vub\ = (3.6 ± 0.4 ± 0.2±lg j) X 10-3 [148]. 

BaBar [149] uses a sample of 206 fb"^ to obtain B{B Tr+tvi) = (1.46 ± 0.07 ± 0.08) x 10^^. 
Recently, the availability of very large data sets at Belle and BaBar have made possible tagged 
analyses, where semileptonic decays are studied in samples where the other B is fully 
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reconstructed, thus defining the event kinematics even more precisely. The first 
implementation relies on the the partial reconstruction of exclusive B — > D^*^£'^U£ decay to tag 
the presence of a BB event [150,151]. 

Belle uses fully reconstructed hadronic tags, achieving the best kinematic constraints and thus 
the highest background suppression capabilities and the most precise determination of q'^. 
Branching fractions obtained with this technique have a bigger statistical error because of the 
penalty introduced by the tag requirement, but the overall error is already comparable with 
the other methods. Table [5] summarizes the present status of the experimental information on 
this decay. Both the total branching fraction for B niui and the partial branching fraction 
for > 16 GeV^ are shown. The latter partial branching fraction is useful to extract \Vub\ 
using unquenched lattice calculations, as this is the only interval where their calculation is 
reliable. 

Table 8: Partial and total branching fractions, in units of lO"'^, for the decay B — > 7r£~^ui. 
Whenever possible, B^ and B~^ data are combined. 



Experiment total 



CLEO [148] 
BaBar [149] 
Belle [150] 
Bellet [152] 
Bellet [152] 
BaBar [151] 



q 



^ > 16 GcV^ 



Method 



combined 



1.37±0.15±0.11 
1.46 ±0.07 ±0.08 
1.38 ±0.19 ±0.14 ±0.03 
1.49 ±0.26 ±0.06 
1.53 ±0.20 ±0.06 
1.33±0.17±0.11 



0.41 ±0.08 ±0.04 
0.38 ±0.04 ±0.04 
0.36 ±0.10 ±0.04 ±0.01 
0.31 ±0.12 ±0.01 
0.39 ±0.12 ±0.02 
0.46 ±0.10 ±0.06 



Untagged Analysis 
Untagged Analysis 
B D^*Hu Tag 

TT^£^i/£ (hadron tags) 
Combined and B~^ tags 
Combined and _B+ tags 



1.34 ±0.06 ±0.05 



0.37 ±0.03 ±0.02 



HFAG ICHEP08 



f - Preliminary results. 



In order to interpret these results, we need theoretical predictions for the form factor /^.(g^). 
This problem can be split into two parts: the determination of the form factor normalization, 
/+(0), and the functional form of the dependence. Form factor predictions have been 
produced with quark models [153] and QCD sum rule calculations [154]. Lattice calculations 
provide evaluations of f+{q'^) at specific values of or, equivalently, pion momenta (p-k)- 
Authors then fit these data points with a variety of shapes. Typically a dominant pole shape 
has been used in the literature. Nowadays more complex functional forms are preferred. 
Becirevic and Kaidalov (BK) [155] suggest using 

f / 2m cb(1 -g) , . 

(1 - gVA/|*)(l - aq-^/M^.) 

where CBM^^, is the residue of the form factor at q"^ = M^*, and Mgi,/a is the squared mass of 
an effective 1~ B*' excited state. Ball and Zwicky [154] propose 

f+{(f')\BZ = - 2nj2 +1 TnTT' *^^^) 

1 — q'^ /M^^, 1 — aq'^/M^^ 
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where the parameters ri, r2, and a are fitted from available data. Lastly, parameterizations 
that allow the application of constraints derived from soft-collinear effective theory (SCET), 
and dispersion relations have been proposed by Boyd, Grinstein, and Lebedev [92], and later 
pursued also by Hill [156]. They define 

where t is {ps — Pw)'^, defined beyond the physical region, t± = {niB ± rriT^)'^, and to is an 
expansion point. BGD use to = 0.65t_. The parameter ax allows the modeling of different 
functional forms, and the variable 

^(t,to) = ^^S ~^ /+~^° (70) 

maps t+ < t < oo onto |z| = 1 and — oo < t < t+ onto the z interval [-1,1]. 

All are refinements of the old ansatz of a simple pole shape, now rarely used. The first lattice 
calculations were carried out with the quenched approximation that ignores vacuum 
polarization effects [157,158]. In 2004, preliminary unquenched results were presented by the 
Fermilab/MILC [159] and HPQCD [160] collaborations. These calculations employed the 
MILC collaboration Nf = 2 + 1 unquenched configurations, which attain the most realistic 
values of the quark masses so far. Using these calculations and the most recent value of the 
partial branching fraction B{B nlu) for > 16 GeV^, shown in Table [H we obtain 
\Vub\ = (3.51 ± 0.096exp ± 0A9th) x 10"^ ^j^i^ ^j^g HPQCD normalization, and 
]T4j,] = (3.70 ± O.Wexp ± 0.37th) X 10"^ with the Fermilab/MILC normalization. Fits to 
experimental data, combining lattice predictions, and dispersion relations reduce the 
theoretical errors. For example, [161] obtains \Vub\ = (3.5 it 0.17 it 0.44) x 10~^, and, more 
recently, [162] obtains \Vub\ = (3.47 it 0.29) x 10~^. The most recent lattice calculation [163] 
performs a simultaneous fit of improved lattice numerical Monte Carlo results and the 12-bin 
BaBar experimental data on \ Vub\f+{q'^) [164] and derives \Vub\ = (3.38 it 0.35) x lO^'^. The ~ 
10% error includes theoretical and experimental errors, not easily separable. 



4.2 Determinations Based on Inclusive B Semileptonic Decays 

Inclusive determinations of \ Vub\ rely on the heavy quark expansion (HQE), which combines 
perturbative QCD with an expansion in terms of l/nib, which accounts for non-perturbative 
effects. Although the possible breaking of the assumption of local quark hadron duality may 
produce unquantified errors, other non-perturbative uncertainties can be evaluated with 
systematic improvements, and their uncertainties are easier to assess than the ones of 
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unquenched lattice QCD or QCD sum-rules. This statement applies to the total charmless 
semileptonic width T^^. In the OPE approach, discussed in the \ Vcb\ section, the observation 
by Chay, Georgi and Grinstein that there are no non-perturbative corrections of order 
^QCD / rrib [93] has inspired the hope that this approach would lead to a more precise 
determination of this important parameter. Before discussing the methods used to relate data 
with theory, it is useful to observe one of the possible expressions available in the literature for 
the total semileptonic width [165], including an exact two-loop expression for the perturbative 
expansion and second-order power corrections. They obtain 



CplVubl'^ml 
r« = I92vr3 + 



-0.768 2.122— 

rrib 



(71) 



-a|(/i) 
3(/u^-A2) 



2 2 

2.158 + 1.019 In — + (1.249 + 2.184 In — + 0.386^/i^/mb) + 0.81 A 



ml 



+ ...} 



where fj, ~ rub is the scale at which as needs to be evaluated, while the scale /x* applies to the 
non-perturbative expansion parameters, namely the b quark mass nib, the chromo-magnetic 
operator, A2, and the kinetic operator, /i^r. It is clear that, even restricting our attention to the 
total width, a precise knowledge of the b quark mass is critical, and considerable theoretical 
effort has been devoted to a reliable extraction of this parameter from experimental 
observables. Other uncertainties, such as the effects of weak annihilation or violations of 
quark-hadron duality, will be discussed later. 

Charmless B meson semileptonic decays constitute only about 1% of the total semileptonic 
width. Thus the big challenge for experimentalists is to identify techniques to suppress this 
large background. For example, the first evidence for B meson charmless semileptonic decays 
came from the study of the end point of the lepton spectrum, where leptons from b — > civ 
processes are forbidden due to the larger mass of the hadronic system formed by the c 
quark [166]. While this was very important first evidence that \Vub\ / 0, very quickly several 
authors pointed out that this region of phase space is ill suited to a precise determination of 
\Vub\ because near the end point the OPE does not converge and an infinite series of 
contributions needs to be taken into account [167-169]. Thus, a large effort has gone into 
developing experimental techniques that would feature a low lepton energy Eg cut and an 
acceptable signal to background ratio. Table [9] summarizes the present status of the \ Vub\ 
determination with this approach. 

Next a whole host of papers proposed alternative "model independent" approaches to measure 
\Vub\ from inclusive decays [165, 173, 174], with the common goal of identifying a region of 
phase space where experimentalists can suppress the 6 — > c background, and where the OPE 
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Table 9: Summary of the experimental \Vub\ determinations using the lepton endpoint; \Vub\ is 
extracted with the "shape function" method [165]. 

Experiment Lepton energy range Cmt (fb~^) AB{Xuli^e) (xlO^) \Vub\ (xlO"^) 

T(45) cont. 

CLEO [170] 2.6 GeV > > 2.1 GeV 9.1 43 3.28 ± 0.23 ± 0.75 3.94 ± 0.46l^J 
Belle [171] 2.6 GcV > > 1.9 GeV 27 8.8 5.72 ± 0.41 ± 0.65 4.74 ± 0.44+gJ 
BaBar [172] 2.6 GeV > > 2.0 GeV 80 9.5 5.72 ± 0.41 ± 0.65 4.29 ± 024±g;| 



works. The first proposal by Bigi, Dikeman, and Uraltsev proposed considering semileptonic 
events where Mx < 1.5 GeV [173]. However, Bauer, Ligeti, and Luke pointed out [175] that 
the kinematic limit nix ~ mf) has the same properties of the lepton end point, and spoils the 
convergence of the OPE; the same authors proposed using dV/dq^ up to = and argue 
that this distribution is better behaved in the kinematic region of interest. This is the 
theoretical foundation of the so called "improved end point" method, which encompasses the 
simultaneous study of Ec and g^. BaBar [176] used this technique, inferring the v 
4-momentum from the missing momentum (pmiss) in the event. There results give \ Vub\ central 
values between (3.88-4.93) x 10^"^ , depending on the explicit model and with total errors of 
RilO% on each value. 

An alternative approach [165], incorporates hadronic structure functions to model the region of 
large hadronic energy and small invariant mass, not well modeled by the OPE, and applies the 
OPE to the kinematic region where the hadronic kinematic variables scale with Mb, and 
smoothly interpolates between them. This approach is commonly referred to as "shape 
function" method, and uses the 7 spectrum in inclusive B — > to reduce the theoretical 
uncertainties. It combines the experimental data on high momentum leptons from B decays 
with the constraints from inclusive radiative decays, to produce a precise value of \Vub\ 
inclusive. 

The study of charmless inclusive semileptonic decays benefits from the use of hadronic tags. 
Belle used the information of the tag momentum to boost the electron into the B meson rest 
frame and to select a sample of high purity. They then reconstruct the i' energy and 
momentum from the measured 4-momentum vectors of the T(45), B tag, lepton, and the 
additional tracks not used to form the tag or the lepton. They then evaluate the invariant 
mass Mx and the quantity P+ = \Ex - Px\, where Px = Pbeam - Ptag -Pl-Pv Both Mx 
and P"*" are smaller for 6 — >^ u transitions. They define different signal region for Pi> 1 GeV: 
P+ < 0.66 GcV/c, Mx < 1.7 GeV/c^, and Mx < 1.7 GeV/c^ combined with q"^ > 8 GeY'^/c^. 
They evaluate the partial branching fractions in each of them, and extract \Vub\ directly from 
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the partial branching fractions, normahzed by corresponding theoretical scale factors i?(A$), 
evaluated with the shape function method [177]. The authors observe that different kinematic 
cuts give different values of \Vub\ and speculate that this may be due to additional theoretical 
uncertainties not completely accounted for. Recently, BaBar has also used the same technique 
in slightly different kinematic regions: P+ < 0.66 GeV/c, Mx < 1.55 GeV/c^, and Mx < 1.7 
GeV/c^ combined with q"^ > 8 GeV^/c^. They use both the shape function method [177] and 
the dressed gluon exponentation method [178]. The results obtained with this approach are 
summarized in Table [TOl 



Table 10: Inclusive \Vub\ determinations with tagged samples. 

Experiment Cint{^~^) |Kfe| x 10^ 

T(4g) cont. 

Belle [171] 253 18 4.09 ± 0.44+^;i]|5 
BaBar [172] 347.4 4.21 ± 0.20l[5:^7 



In order to provide predictions that are most suited to different experimental cuts, theorists 
have made available the triple differential width ^^^^Ji^g^ given by 

G||K6|'_^^2^^_p£;|_2goS,+^]^2+g'(2£;£-(?o)W^3}x^(go-^£)^)^(^^)^(g')^(go- 



(72) 

where qo is the energy of the lepton-zv pair and Ei is the energy of the charged lepton in the B 
meson rest frame, and VF1-3 are the three structure functions relevant if we assume massless 
leptons. Reference [179] computes the functions Wi{qo,q'^) as a convolution at fixed q'^ between 
non-perturbative distributions Fi{Kj^,q^] ji) and the perturbative functions Wf^"* {qo,q^) 



wMo^q^) = J dk+F,ik+,q'■,^,)w^' [go - ^ (1 



,2 



mbMB 



■,mb] . (73) 



Perturbative corrections to the structure functions VFi_3 are now known up to order 
C(a2/3o) [180] and power corrections are included through O(m^). The separation between 
perturbative and non-perturbative physics is set by a cut-off scale /U = 1 GeV, in the "kinetic 
scheme" [142], which takes input parameters from fits to the b ^ c moments. Alternative 
choices of kinematical variable have been used, for example, in the SCET approach, the 
variables = \Ex — Px\ and P~ = \Ex + Px\ are used. In particular, the "shape function" 
method [165] provides theoretical expression for the triple differential d^T / dEidP~^ dP" and 
relates it to moments of the shape function extracted from B ^57. Finally, the approach 
originally proposed by Gardi [181] uses resummed perturbation theory in momentum space to 
provide a perturbative calculation of the on-shell decay spectrum in the entire phase space. 
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Figure 25: Summary of Inclusive \Vub\ determinations. 

The metliod used, dressed gluon exponentiation, is a general resummation formalism for 
inclusive distributions near a kinematic threshold. Finally, another model based on soft-gluon 
resummation and an analytical time-like QCD coupling has been proposed [182]. Clearly a lot 
of work has gone into bringing to fruition the original promise that the inclusive \Vub\ 
determination is a more precise method to determine this important CKM parameter. 
However, Figure [25] shows that different methods provide central values of \Vub\ that often 
differ beyond the stated errors. Moreover, these estimates are generally significantly higher 
than the value of \ Vub\ extracted from B ■nivf^ and the one that is obtained from global 
unitarity triangle fits. 

An effect that can influence the inclusive \Vub\ is the contribution due to topologies where the 
incoming hq pair annihilates into a W boson, with a concomitant soft gluon emission. These 
contributions are known as "weak annihilation" and they appear as a delta function at the end 
point of the lepton spectrum. This effect introduces a difference between and i?^ 
semileptonic decays. 

A first attempt to quantify the annihilation effect was performed by CLEO [183], studying the 
high contributions to charmless semileptonic decays. They used inclusive data on a sample 
including both charged and neutral B mesons. They used a variety of models to set the limit 
on the fractional contribution of annihilation diagrams, FvyA/Ffe^+u < 7.4 % at 90 % confidence 
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level. BaBar has presented a measurement of B{B^ X^ii/i) at the end point region which, 
combined with previous inclusive measurement, sets the limit 

^^^/^'^-^ < /^^(2.3 - 2.6) ^^^^ 

at 90 % confidence level, where fwA{'2-3 — 2.6) represents the fraction of weak annihilation 
rate contributing in the momentum interval Ap£ comprised between 2.3 and 2.6 GeV used in 
this analysis, believed to be close to unity. 

Finally, Voloshin [184] has suggested that the difference in semileptonic widths of the and 
Dg mesons can assess non-factorizable terms through the relationship 



Toic^slu) ' V 0.22 GeV 

where Tq{c siu) is the bare parton semileptonic rate c — > sii', fo is the D decay constant, 
and S"'^ — i?2 is the non-factorizable term that affects charmless semileptonic decays as well. 
Using this formula with the CLEO-c results for inclusive semileptonic decay [185], the sum 
of Dg exclusive semileptonic decays [186], and mdmc) = 1.268 GeV [187], we obtain an 
estimate of non-factorizable effects equal to (9zb5)%. Note that annihilation effect induced 
uncertainties are estimated to contribute 1.6% in the BLNP approach, 1.5% in the DGE 
approach, and (-1-0.0,-3.1)% in the GGOU approach. 

In summary, the values if \ Vub\ extracted with the inclusive method have quoted errors between 
5.5% and 10%, with central values that change well outside these uncertainties and are all 
considerably larger than the \ Vub\ value obtained from the exclusive branching fraction 



5 Rare B Decays 

In general, we define as "rare" B decays processes that are suppressed at tree level. They are 
interesting because they are typically mediated by loop diagrams which may be characterized 
by a matrix clement whose strength is comparable to components with similar Feynman 
diagrams, where new particles appear in loops. Thus evidence for new physics may appear 
either through enhancements in branching fractions relative to the Standard Model 
expectation, or through interference effects. 
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Table 11: Summary of recent theoretical evaluations of the decay constant fs- 



Method /b (MeV) 

HPQCD (unquenched) [188] 216 ± 9 ± 19 ± 4 ± 6 
LQCD (chiral extrapolation) [189] 209.4 ± 9.7 ± 1.0 
FNAL-MILC-HPQCD [190] 197 ifc 6 ifc 12 



Table 12: Summary of experimental determinations of B{B tu) 

Experiment Method B{B tu) x 10^ 

Belle [191] Hadronic tag sample 1.79+^:^^+^:|'j' 

Belle [192] Semileptonic tag sample 1.65t[!|^l[!|^ 

BaBar [193] Hadronic tag sample l-8lo:8 ± 0.4 ± 0.2 

BaBar [194] Semileptonic tag sample 1.8 ± 0.8 ± 0.1 

Average 1.73 ±0.35 



5.1 B^Tv 

The decay B ^ tv \s affected by two quantities of great interest, the quark mixing parameter 
\Vub\ and the pseudoscalar decay constant Jb- In fact, the leptonic branching fraction is given 
by 

B{B^ ^ r^u) = j^i _ /||K,|V^, (76) 

where uij- is the mass of the r lepton, and tb is the charged B lifetime. Theoretical predictions 
for Jb are summarized in Table [TTJ The most recent value exploits the full machinery of 
unquenched Lattice QCD, and has a precision of 7%. 

Belle performed the first measurement of the branching fraction B{B^ t^v) using a tagged 
sample of fully reconstructed B decays [191]. They later used a tagged sample of semileptonic 
decays [192]. BaBar also published studies performed using fully reconstructed B decays [193] 
and semileptonic tags [194] . The results are listed in Table [T2l The resulting average 
experimental value B{B tv) is (1.73 ± 0.35) x 10""^. 

The Standard Model prediction for B{B tv) is (O.77+0I2) x 10"^, obtained using the 
exclusive value of \Vub\ [163] and the most recent Lattice QCD value of /_b [190]. The range 
has been obtained by calculating the Standard Model branching fraction using the values 
iKbl = WuhUxd ± cr(lKfel) and Jb = fB± cr{fB). An intriguing possibility is the enhancement 
of this branching fraction due to new physics. For example in a 2-Higgs-doublet model 
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Figure 26: B{B — > tv) as a function of the effective scalar coupling gg. The shaded region 
is the measure data ztlo" The range of theoretical predictions is obtained by assuming \ Vub\ = 
(3.38±0.35) X 10-3 [163], and fs = (0.197±0.06±0.012) MeV [190]. (The curves are theoretical 
and explained in the text. 



(2HDM) [195] based on supersymmetric extensions of the Standard Model, the ratio 



where Qs is the effective scalar coupling. Fig. [26] shows the predicted value of B{B tv) as a 
function of Qs and the band representing the measured value, which constrains Qs to be less 
than 0.1 or between 1.8 and 3.9. A better understanding of the \Vuh\ and Jb inputs is 
necessary to improve on this estimate. 



experimental value for this parameter is (41.6 it 11.7 it 5.2)% [7], which constrains gs to be 
< 1.79. Thus these two measurements provide complementary constraints, which have very 
different systematic uncertainties. Thus an improvement in the knowledge of leptonic and 
semileptonic branching fractions involving a r lepton in the final state is very important. 

5.2 Radiative B Decays 

Radiative B decays involving h — > s{d)^ transitions are very sensitive to new physics processes. 
These processes are ideals for indirect searches for physics beyond the Standard Model, such as 
two-Higgs doublet models , supersymmetric theories and extended technicolor scenarios [196]. 
Hence, comparison of results from these theories with experimental measurements places 
constraints upon new physics. Moreover, h — > Xg-^ is an ideal laboratory for studying the 




(77) 



Note that a quantity that is sensitive to the same 



coupling is the ratio -wr^zrfjr^- present 
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dynamics of the 6-quark inside the B meson: since the motion of the 6-quark inside the B 
meson is universal, information gained from a measurement of the energy spectrum of the 
emitted photon in this transition is apphcable to other processes, for instance semileptonic 
decays. 

In general, in the OPE, the amplitude A for a process can be expressed as sum 

A = (Hefi) = E Mw){Q^{^^)), (78) 

i 

where the Qi are local operators, the Cj are Wilson coefficients, and fi is the QCD 
renormalization scale [197]. The B — > branching fraction directly probes the Wilson 
coefficient C7. However, some new physics contributions may leave the B — > rate 
unaltered, with changes to the C7 amplitude. In this case the direct CP asymmetry is sensitive 
to new phases that may appear in the decay loop [198] and the B Xs£^i~ transition may 
provide information on the sign of the amplitude, since it also probes Cg and Ciq. 

The B Xs'j branching fraction, as a function of a photon energy cut-off £^0) is related to the 
B Xcvi transition by 



B{B^ Xsj ; E^> Eo GeV) = B{B ^ X^eu] 



exp 



Gap 



ttC 



X [P{Eo) - NiEo)] (79) 



Vr 



cb 



V 



ub 



TjB ^ X^ef?) 
T(B ^ X^eu) 



where the perturbative corrections P{Eq) are defined as 



VtlVd 



V 



cb 



6ae_ 



vr 



-P{Eo 



Vcb 


2 


T(6 S-f)E>Eo 


Vub 




T{b cev) 



(80) 



LO 



(1 + 0{as)NLO + 0{a)NLO + O 



and the non-perturbative N{Eq) terms are 0{^)^q,^nlo+0{-^)];^q,^nlo 

^)NLO+m^° ■ Note that the minimum value of £'0 for which this relation is valid is 1.6 
GeV. The theoretical error from NLO perturbative calculations is about 10%, dominated by 
the renormalization scale dependence and a charm quark mass uncertainty. This uncertainty is 
mainly due to the change from the ratio in the pole scheme m^°'*^/m^°'^ = 0.29 it 0.02 to the 
ratio m^^/m^^ = 0.22 it 0.04 in the MS scheme (see section [3.0. ip . The large ruc dependence 
is due to the fact that nic first enters at NLO matrix elements. Hence, the natural scale at 
which nic should be renormalized can only be determined by dedicated calculations at NNLO. 
Two predictions for the branching fraction at NNLO have been given: 

B{B ^ X,7 ; £^>i.6 Gev) = (3.15 ± 0.23) x 10-^ (81) 

= (2.98 ± 0.26) X 10"*^, 
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from [199] and [200], respectively. Both agree with the world average in Table [T3l The 
theoretical error of 7% is obtained by quadratically adding the uncertainty due to 
non-perturbative corrections (5%) [201], the uncertainty due to missing higher-order 
corrections (3%), the ttt-c interpolation ambiguity (3%) and the parametric uncertainty (3%). 
Not included are some known NNLO and non-perturbative corrections. The size of these 
neglected contributions is about 1.6%, which is smaller than the present theoretical error. 
While progress in the calculation of the perturbative corrections is expected in the future, the 
uncertainty on the non-perturbative effects will not be easily reduced as they are very difficult 
to estimate. 

Experimentally, two methods are used to extract the b — > Xg^ signal: the fully inclusive 
method and the semi-inclusive method. In the fully inclusive method, events containing a hard 
photon consistent with B Xg'j are selected. In this method the subtraction of a very large 
background, primarily from qq continuum events, is the main issue. There are several options 
to reduce this large background. Particularly effective in the suppression of continuum 
backgrounds is the requirement of a high momentum lepton, tagging the semileptonic decay of 
the accompanying B meson. Alternatively, a fully reconstructed B can be used to identify B 
decays. In this measurement the large background increases towards the lower photon energy, 
making it impossible to measure the full photon spectrum. Hence, all measurements require a 
minimum photon energy where the signal to noise ratio is still acceptable. The total branching 
fraction is then obtained by extrapolating the signal to the full phase space. This 
extrapolation is based on theoretical models and is an irreducible source of systematic error. It 
has been argued that the energy scale A = mf, — 2£'™™ is significantly smaller than mi, and 
therefore the above extrapolation error is underestimated. Therefore it is important to 
measure the b ^ photon energy spectrum as precisely as possible. So far we have been able 
to lower the minimum photon energy to 1.7 GeV. 

In the semi- inclusive method, the B Xg'j branching fraction is determined by summing up 
exclusive modes, with an extrapolation procedure to account for the unobserved modes. The 
semi-inclusive analysis provides a more precise photon energy in the i?-meson rest frame and 
usually provides the flavor of the decay (i.e. 6 or b). The extrapolation to account for the 
unobserved modes is is the key issue in this approach; it is based on isospin symmetry and the 
S-meson Monte Carlo hadronization model, which is hardly precise and reliable. All the 
available measurements are summarized in Table [131 Note that the CLEO measurement on the 
branching fraction includes b ^ dj events, which is expected to be quite small 
~ iVtd/Vtsl^ « 5%. 

At the parton level in the two-body decay b — > 57, the energy of the photon is w nih/2 in 
the 6-quark rest frame. However, the b Xg'j is not mono-chromatic, due to several effects 



58 



including the width of the Xg mass distribution, gluon emission and the Fermi motion of the 
6-quark in the B meson. The non-perturbative Fermi motion effects will raise the photon 
energy above mb/2, while the gluon emission will give a long low-energy tail. The low-energy 
tail can be described by HQE a part non-perturbative effects that can be modeled by 
process-independent shape functions described by a few parameters, e.g. the h quark mass and 
the Fermi momentum (/i^) that are considered to be universal in the following inclusive B 
decays: B X^'j, B Xgii, B X^lf and B Xuiv. Recently a significant effort was 
made to combine all available data for B ^ Xgj and B Xciv to determine the HQE 
parameters. Using these parameters, B ^ Xs"f branching fraction results are combined 
together to provide a rather precise branching fraction. The method to combine these 
measurement has been provided by Limosani [202] and it is different from the method used by 
HFAG, as explained below. 

The extrapolated branching fractions do not include the published model uncertainties and the 
uncertainty on the extrapolation factor. These errors are included in the average by 
recalculation in the framework of a particular ansatz. Correlations between different 
measurements have been ignored. The parametric error on rrn, is evaluated by varying m;, 
within its uncertainty. The world average has been calculated taking into account the 
correlations, when available between partial branching fraction measured at different photon 
energy thresholds for each single analysis. In the branching fraction average we use the energy 
threshold of each photon energy spectrum that corresponds to the optimal overall uncertainty 
on the full rate after extrapolation. This differs from the HFAG method, which uses the lowest 
energy thresholds for each spectrum measurement. The HFAG method penalizes analyses that 
quote measurements at low Emin, which suffer from larger systematic uncertainties. 

Two different calculations were used to extrapolate the measured partial branching fractions 
down to a photon energy lower threshold of 1.6 GeV. The extrapolation factors were 
determined using m^^ = 4.243 ± 0.042 GeV as input [178], or with mf^ = 4.63 ± 0.04 GeV 
and nl = 0.272l|j;[|^^GeV2 as input [200]. The two world averages are listed in Table US and 
they agree within their theoretical errors. 

The agreement between the measured B — > Xs'j branching fraction and the theoretical 
prediction constrains various new physics scenarios. One of the most popular examples is the 
lower bound on the type-II charged Higgs mass, since it always constructively interferes with 
the SM amplitude. The current limit is around 200 GeV for any tan /3 if no other destructive 
new physics contribution exists [203]. This limit is significantly higher than the direct search 
limit. 

^For the time being only Belle provides these correlations. 
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Table 13: Measured branching fractions, minimum photon energy, and branching fractions for 
Emin = 1-6 GeV photon energy for b sj. The third error is the model uncertainty quoted 
by the experiment. Two world averages are calculated extrapolating all the branching fractions 
using two different theoretical calculations [178] and [200]. The Belle semi- inclusive measurement 
sums up 16 modes, BaBar 38 modes. The Belle inclusive analysis calculates the branching 
fraction for different photon energy cut-offs. 



Experiment 


-'-'mm 


B (10-«) 


SModif.(10-«) [200] 


i3Modif.(10-^ 


) [178] 






at Ejjiiii 


S^in = 1.6 GeV 


Emin = 1-6 


GeV 


CLEO(Incl.) [204] 


2.0 


306 ± 41 ± 26 


337 


323 




Belle(Semi-inc) [205] 


2.24 


336 ± 53 ± 42 


496 


434 




BaBar(Semi-inc) [206] 


1.9 


327 ± 181:11 


354 


337 




BaBar (Incl.) [207] 


1.9 


367 ± 29 ± 34 


397 


378 




BaBar(recoil) [208] 


1.9 


366 ± 85 ± 60 


396 


377 




Belle(Incl.) [209] 


1.7 


331 ±19 ±37 


337 


333 




Belle(Incl.) [209] 


1.8 


324 ± 17 ± 24 


339 


329 




Belle(Incl.) [209] 


1.9 


312 ± 15 ± 16 


338 


321 




Belle(Incl.) [209] 


2.0 


294 ± 14 ± 12 


334 


310 




Average 


1.6 


350 ± Ue^p ± 5™, ± 8^2 


322 ± Uexp 


± ^mj, 








xVndf = 7.3/8 


X^/ndf = I 


i3/8 



New physics contributions may leave the B — > X^^ rate unaltered, with new phases appearing 
in the decay loop [198]. Since the Standard Model prediction of the CP asymmetry {Acps^y) is 
zero in the limit of U-spin symmetry [210], significant non-zero values would be evidence for 
new phenomena. AcPs'y has been measured by Belle, BaBar and CLEO. The values of all such 
measurements are listed in Table [TH All measurements are in agreement with the Standard 
Model predictions. BaBar also report a CP-asymmetry for b ^ [d + s)j of 
AcPs'y = —0.11 ± 0.12 ± 0.02 [207], measured using an inclusive analysis with a lepton tag. 

The B-factories also attempted to measure the polarization of the photon in the 6 — > 57 
transition, which can provide a test of the Standard Model, which predicts the photon to be 
mainly left-handed [198]. The measurements rely on either the exploitation of the — B^ 
interference [198] or the interference between higher-Kaon resonances decaying into 
K-KTT^ [211]. However, both Belle and BaBar [212] measurements are not precise enough to pin 
down the photon polarization. The list of established decay modes of the type B K^*^Xj, 
where X is one or more flavorless mesons, is increasing. These channels are interesting as they 
can provide a measurement of the polarization of the photon. Many decay involving kaon 
resonances modes (e.g. i^*(892), (1270), ii'|(1430)) are now well established and their 
branching ratios agree with the Standard Model predictions. 
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Table 14: CP Asymmetries for exclusive and inclusive 6 ^ 57 transitions for , and B^ / B^ 
admixture. 



Experiment 




Acp 


BaBar [213] 




-0.013 ±0.036 ±0.010 


Belle [214] 




-0.015 ±0.044 ±0.012 


CLEO [215] 




0.08 ±0.13 ±0.03 


Average [99] 




-0.010 ±0.028 


Experiment 




Acp 


BaBar [216] 




-0.011 ±0.030 ±0.014 


Belle [217] 




0.002 ± 0.050 ± 0.030 


CLEO [218] 




-0.079 ±0.108 ±0.022 


Average [99] 




-0.012 ±0.028 



Table 15: Summary of upper limits for B ^ ii decays. 



Mode 






Experiment UL 9' 


0% CL (xlO«) 


B^^ 






BaBar [219] 


5.2 


B^^ 






CDF [220] 


1.5 


Bs^ 






CDF [220] 


4.7 


Bs^ 






DO [221] 


7.5 


B^^ 






BaBar [222] 


11.3 


B^^ 






Belle [223] 


19 


B^^ 






CDF [224] 


8.3 


Bs^ 






CDF [224] 


28 



5.3 B^U 

Flavor changing neutral current (FCNC) decays proceed only through loop diagrams in the 
Standard Model, and are further reduced by helicity and Glashow-Iliopoulos-Maiani (GIM) 
suppression. Thus even the largest branching fraction predicted by the Standard Model, 
Bsm{Bs — > ^J'~^fJ'~) ~ 4 X 10~^, is very tiny. New physics can enter either at tree level, for 
example through the presence of an additional Z' boson, or through new particles appearing in 
loops, and can increase the actual value of these branching fractions. Decay modes involving 
two muons in the final state are particularly interesting because they are amenable to 
measurements in experiments operating at hadronic machines. In fact the best upper limits 
come from CDF and DO, which are reaching a sensitivity comparable to the SM prediction. 
Table [15] shows a summary of the present experimental values. 
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Table 16: Summary of inclusive and exclusive b Xgii branching fractions. 



Mode 



BaBar [225] 



Belle [226] 



CDF [227] 



Average [99] 



K+i+r 

K+e+e~ 

K*+i+i- 
K*+e+e- 




0.59 ±0.15 ±0.04 



0.49 ± 0.05 
0.52 ±0.07 




-0.27 



5.4 B ^ XsU 

The physics of i? — > Xgii is governed by the Wilson coefficients Cy, Cg and Cio, which 
describe the strengths of the corresponding short-distance operators in the effective 
Hamiltonian, i.e. the electromagnetic operator O7 and the semileptonic vector and axial-vector 
operators O9 and Oio, respectively [197]. The Wilson coefficients are experimental observables. 
Contributions from new physics appear in the experiment as deviations from the SM values, 
which have been calculated to next-to-next-to-leading order (NNLO). The experimental 
knowledge on the Wilson coefficient C7 comes from the inclusive b — > Xs'y branching fraction, 
which determines its absolute value to about 20% accuracy, but not its sign. The partial 
b XgU decay rate in the lepton invariant mass range below the J/ij: resonance is sensitive to 
the sign of C7. Measurements of the inclusive b Xgii decay rate have been published by 
Belle [228] and BaBar [229] (see Table [T6l who also report a direct CP asymmetry compatible 
with zero, see Table [iTl 

The inclusive B Xgii branching fraction, which constrains Cg and Cio [230], gives no 
information on the individual signs and magnitudes of these coefficients. To further pin down 
the values of these coefficients, it is necessary to exploit interference effects between the 
contributions from different operators. This is possible in B ^ Xgii decays by evaluating the 
differential inclusive decay rate as a function of the lepton invariant mass, = m[£i)'^, or by 
measuring the forward-backward asymmetry in the exclusive decay B K*ii. 

The forward-backward asymmetry in S — > K*U, defined as 



an ideal quantity to disentangle the Wilson coefficients Cj since the numerator of Eq. [82] can 



(82) 



is a function of q'^ and of 9^1-^ the angle between the negative lepton and the B meson. It is 
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Table 17: CP Asymmetries for exclusive and inclusive b — > Xg££ transitions for S^, and 
admixt ure. 



Mode 


BaBar [231,232] 


Belle [226] 


Average [99] 




-0.07 ±0.22 ±0.02 


-0.04 ±0.10 ±0.02 


-0.05 ±0.09 


K*+ee 




-0.13 ±0.17 ±0.02 


-0.13 ±0.17 


K^U 




-0.08 ±0.12 ±0.02 


-0.08 ±0.12 


K*^U 




-0.10 ±0.10 ±0.02 


-0.10 ±0.10 


K*£i 


0.03 ±0.23 ±0.03 




0.03 ±0.23 


sU 


-0.22 ± 0.26 ± 0.02 




-0.22 ± 0.26 


K*U 


0.03 ± 0.23 ± 0.03 




0.03 ± 0.23 



be expressed as 

^ 



Re(C9)Fi ± ^CyFs 



(83) 

L y- J 

where ^ is a function of q^, and F12 are functions of form factors. It is straightforward to 
determine Cio, Ke{Cg) and the sign of C7 from the ApB distribution as a function of q^, using 
the value of ICyj from B — > Xg-f and a few more assumptions: phases of Cio and C7 are 
neglected, and higher order corrections are known. These assumptions should be examined by 
comparing the results with the inclusive B Xgll differential branching fraction as a function 
of gr^, since it is also sensitive to Cg and Cio in a different way. Most of the eight individual 
B —>■ K{K*)U modes have been established. Both experiments have searched for asymmetries 
with respect to lepton flavour. Currently available data clearly favour a negative sign for C7, 
as predicted by the Standard Model. 

In some SUSY scenarios the sign of the 6 ^ 57 amplitude (C7) can be opposite to the SM 
prediction, while the transition rate may be the same. Within the SM there is a zero crossing 
point of the forward-backward asymmetry in the low region, while it disappears with the 
opposite sign C7 if the sign of Re(C9) is the same as in the SM. In another model with SU{2) 
singlet down-type quarks, tree-level Z fiavor-changing-neutral-currents are induced. In this 
case, the larger effect is expected on the axial- vector coupling (Cio) to the dilepton than on 
the vector coupling (Cg). Because the forward-backward asymmetry is proportional to the 
axial-vector coupling, the sign of the asymmetry can be opposite to the SM. The same new 
physics effect is also a possibility for B — > 0i^5 where anomalous mixing-induced CP violation 
can occur. 

Belle attempted to extract Cg and Cio from ApB '^^ B ^ K*U with 357 fb"^ data [233]. The 
Belle analysis constrains the sign of the product CgCio to be negative as in the SM. In this 
study higher order QCD correction terms are assumed to be the same as in the Standard 
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Model. Only the leading order terms, Cg and Cio, are allowed to float in the fit to the data. 

bmce is an all charged particle final state, LHCb may be able to measm'e the 

zero crossing point with a better precision than a B factory. However, a model independent 
analysis requires measurement of the forward-backward asymmetry in S — > Xgii, which is only 
possible at a super-S factory. 

The precision of the B K{K*)U branching fractions is dominated by the theoretical errors, 
which have large model dependent irreducible uncertainties in the form-factors, which can be 
as large as a factor two. These uncertainties are much smaller if the ratio 

= B{B K{K*)fi+fi~)/B{B K{K*)e+e-) is measured. In the SM, the B Ke+e- 
and B — > Kfi'^fi" branching fraction are almost equal, a part from a small phase space 
difference due to the lepton masses. However, the B — > K*e~^e~ branching fraction is expected 
to be larger than the B K*fj,^fj,^ branching fraction, due to a larger interference 
contribution from B K*'j in i? — > K*e^e~ . In new physics models with a different Higgs 
sector than that of the SM, scalar and pseudo-scalar types of interactions may arise in b ^ sU. 
Therefore, Rk = B{B K^^ ^^)/B{B i^e"*"e~) is an observable that is sensitive to new 
physics. For example, a neutral SUSY Higgs contribution can significantly enhance only the 
B — > K{K*)ii^^~ channel at large tan/3. The current world average for the branching 
fractions gives: Rk = 0.96 it 0.26, dominated by the statistical error. 

Finally, the transition h — > suT> is theoretically much cleaner than its charged counterpart, 
thanks to the absence of the photon penguin diagram and hadronic long-distance effects 
(charmonium resonances). From the experimental point of view these modes are very 
challenging due to the presence of the two neutrinos. Searches for such modes at the B 
factories use a fully reconstructed B sample or semileptonic tagged mode. Both Belle and 
BaBar report limits for these modes. 

6 Current Status of Overall CKM fits 

We have seen that there are no startling departures from the Standard Model. Yet, we need to 
formulate quantitatively how restrictive the measurements are. First of all, we can 
demonstrate that the data are consistent among the different measurements by plotting the 
extracted values of p versus rj along, with their experimental and theoretical errors, shown in 
Figure [271 from many of the measurements described above, also including the measurement of 
CP violation in the neutral kaon system, designed as "ei^;." Physics beyond the Standard 
Model could be revealed by seeing a particular measurement that is not in agreement with the 
others. The width of the bands are the one-standard deviation combined overall errors. For 
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many of the measurements, e.g. \Vub\, ^k, Am^ and Am^, the theoretical errors dominate and 
thus how they are treated is important. The CKM fitter [11] group generally treats the errors 
more conservatively than the UT fit group [20]. The former uses a frequentist statistical 
approach, while the later a Bayesian approach. The Bayesian errors are smaller, but if the 
question is one of seeing new physics or not, a conservative approach may be better. The 

1.5 



1.0 
0.5 
1?^ 0.0 
-0.5 



-1.0 



-1.5 
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P 

Figure 27: The projections in the Jh-rj plane from measurements of the angles and sides of the 
unitarity triangle. The colored area at the apex shows the 95% confidence level point. (From 
the CKM fitter group [69]). 

results for the Standard Model CKM parameters from both groups are listed in Table [181 The 
Table 18: Values of Wolfenstein parameters from CKM fits assuming the Standard Model 



Group 


P 




A 


CKM fitter 


n -|or-+U.U33 

U.iOO_g 016 


0.345+ni^ 0.795+H?^ 


0.2252 ± 0.0008 


UT fit 


0.155±0.022 


0.342±0.014 





central values of p and rj are in agreement. The value of A is in agreement with that found by 
Blucher and Marciano [235]. The fact that there is an overlap region common to all the 
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measurements shows no obvious deviation from the Standard Model, which does not predict 
values for p and fj, but merely that the various measurements must find the same values. 

To search for deviations from the SM it is useful to compare measurements involving CP 
violating loop diagrams, where imaginary parts of amplitudes come into play, with 
measurements of the sides of triangle. The UT fit group shows in Figure [28] a comparison in 
the p-rj plane of angle measurements with those of the sides of the unitarity triangle. The 
overlap regions in both cases are consistent. 




CP violating ^ CP conserving ^ 

Figure 28: Constraints at 68% confidence level in the ^-r/ plane. On the left side are mea- 
surements of the CP violating angles a, (3 and 7 and ex in the neutral kaon system. On the 
right side are measurements of |V^;,|/V^6| and the ratio of Bs to mixing. (From the UT fit 
group [234]). 



It remains to set limits on New Physics. We allow for NP to contribute to the mixing 
amplitude of either Bd or Bg mesons, which we consider separately. Then the amplitude of this 
second order (Ai^ = 2) flavor changing interaction can be expressed as 



n 



SM+NP 
AF=2 



Bg) = [Re(A,) + ilm(Ag)] • {Bg 



n 



SM+NP 
AF=2 



Br, 



(84) 



where q specifies the type of neutral B meson, and Ag gives the relative size of NP to SM 
physics. (The SM point has Re(Ag)=l and Im(Ag)=0.) 



To set limits it is necessary to make some simplifying, but very reasonable, assumptions [236]. 
We assume that there is no NP in the tree level observables starting with the magnitudes \ Vud\, 
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iKisli {Vcbli and \ Vub\- Since since 7 is measured using tree level decays we include it in this 
category. Now we incorporate the measurement of /3 by noting that a = vr — 7 — f]'^^^^ allowing 
that the measured = fj^^ — Arg(A(^). Also used are measurements of the semileptonic 

asymmetries (see Eg. [237], and the constraint equation AF^ = cosc/i^Arf^^ [238]. The 
allowed range in the p-rj plane for NP via B mixing is shown in Figure [29l 




-2-10123 -2-10123 

Re Re Ag 



Figure 29: Constraints on the imaginary part versus the real part of a New Physics amplitude 
proceedings via B mixing for (left) B^ mesons and (right) Bg mesons. The SM point is at (1,0). 
The contours shown are all at 68% confidence level. (From the CM fitter group [69]). 

In both case quite large values of the real part of NP amplitudes are allowed. In the Bg case a 
large negative value of the imaginary part is also preferred, but again the effect is only 2.7a. 
We can conclude that we need to do a lot more work to constrain NP in these transitions. 

7 Conclusions 

The goals of our studies are to elucidate the features of the electroweak interactions and the 
fundamental constituents of nature. The b was discovered at fixed target experiment in a 
proton beam, by the first observation of narrow T states in 1977. Experiments at several e~^e~ 
colliders have increased our knowledge enormously. CLEO provided the first fully 
reconstructed B mesons and measurements of their mass. The Mark II and MAC experiments 
at PEP gave the first measurements of the b lifetime, which were confirmed by PETRA 
experiments. ARGUS first measured B^ mixing. CLEO, ARGUS and the LEP experiments 
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competed on the first precision measurements of the CKM parameters Vcb and Vub- The two 
most recent, BaBar and Belle, have provided measurements of CP violation. 

The hadron collider experiments have measured Bg mixing and are starting to probe CP 
violation in the Bg system. The torch will soon be passed to first hadron collider experiment 
designed to measure b decays, LHCb. The ATLAS and CMS experiments will also contribute. 

In the 31 years since the discovery of the b quark we have learned a great deal about the 
physics of 6-flavored hadrons. We have seen that the Standard Model describes most decays 
well, yet studies continue as any new particles that we believe must exist have influence, 
especially on rare and CP violating B decays. We expect that future experiments will learn a 
great deal about these new particles either by observing their effects on B decays, or by seeing 
that they do not have observable effects. The latter case is one that experimentalists do not 
like, but it can be just as important in firming up our understanding of the basic constituents 
of matter and the forces that operate between them. 

Acknowledgements 

We thank Phillip Urquijo and Antonio Limosani for useful discussions concerning the CKM 
elements. M. Artuso and S. Stone thank the U. S. National Science Foundation for support. 

References 

1. N. Cabibbo, Phys. Rev. Lett. 10, 531 (1963). 

2. J. H. Christenson, J. W. Cronin, V. L. Fitch and R. Turlay, Phys. Rev. Lett., 13, 138 
(1964). In CP, C stands for charge conservation invariance and P for parity invariance. 
The operation of C takes a particle to anti-particle state and the operation of P changes 
left to right. See the article by John Ellis in the CERN courier 
http://cerncourier.com/cws/article/cern/28092, Also see Yosi Nir in 
http://www.symmetrymagazine.org/cms/?pid=1000194, Streaming video of several talks 
is available at 

://www- project. slac.stanford.edu/streaming-media/CPViolation/CPViolation.htm 

3. An excellent resource is the book by I. I. Bigi and A. I. Sanda, "CP Violation," 
Cambridge Monographs on Particle Physics, Nuclear Physics and Cosmology (No. 9), 
(1999). 



68 



4. J. J. Aubert et al, Phys. Rev. Lett. 33, 1404 (1974); J. E. Augustin et al, Phys. Rev. 
Lett. 33, 1406 (1974). 

5. J. D. Bjorken and S. L. Glashow, Phys. Lett. 11, 255 (1964); S. L. Glashow, J. Iliopoulos 
and L. Maiani, Phys. Rev. D2, 1412 (1972); S. L. Glashow, "Charm: An Invention 
Awaits Discovery," IV Int. Conf. on Experimental Meson Spectroscopy, Apr, 1974, 
Boston, Mass. Published in AIP Conf. Proc. 21, 387 (1974). 

6. M. Kobayashi and T. Maskawa, Progress in Theoretical Physics 49, 652 (1973). 

7. Particle Data Group, C. Amsler et al, Phys. Lett. B 667, 1 (2008). 

8. S. L. Glashow, Nucl. Phys. 22 579, (1961); S. Weinberg, Phys. Rev. Lett. 19 1264, 
(1967); A. Salam in Elementary Particle Physics (Nobel Symp. N.8), Ed. N.Svartholm, 
Almquist and Wiksells, Stockholm (1968), p. 367. 

9. To view a cartoon of the mechanism of Higgs mass generation see 
]http://www.exploratorium.edu/origins/cern/ideas/higgs.html 



10. L. Wolfenstein, Phys. Rev. Lett. 51, 1945 (1983). 

11. J. Charles et al. Eur. Phys. J. C41, 1 (2005); see also 



http : // www. slac . Stanford . edu / xor g / ckmfitt er / ckm_intro. html 



12. R. Aleksan et al, Phys. Rev. Lett. 73, (1994) 18; J. P. Silva and L. Wolfenstein, Phys. 
Rev. D55, 5331 (1997). 

13. M. A. Srednicki, "Particle Physics and Cosmology: Dark Matter," Amsterdam, 
Netherlands, North-Holland (1990). 

14. D. Hunter, G. D. Starkman and M. Trodden, Phys. Rev. D66, 043511 (2002). 

15. J. Hewett and M. Spiropulu, Ann. Rev. Nucl. Part. Sci. 52, 397 (2002). 

16. A. Sakharov, JETP, 5, 32 (1967). 

17. M. B. Gavela, P Hernandez, J. Orloff and O. Pene, Modern Phys. Lett. A 9, 795 (1994), 
and references contained therein. 



18. http:/ /www.esi- topics.com/brane/interviews/DrLisaRandall.html 



19. http://www.slac.stanford.edu/econf/C040802/lec_notes/Lykken/default.htm 



20. See http://www.utfit.org/ 



21. H. Albrecht et al. (ARGUS), Phys. Lett. B 192, 245 (1987). 

69 



22. C. Albajar et al. (UAl), Phys. Lett. B 186, 247 (1987), erratum-ibid B 197 565, (1987). 



23. R. Akers et al. (OPAL), Z. Phys. C 66, 555 (1995). 

24. E. Barberio et al. (HFAG), "Averages of b-hadron and c-hadron Properties at the End of 
2007," larXiv:0808.1297. [hep-ex]. 

25. M. Gaillard and B. Lee, Phys. Rev. D 10, 897 (1974). 

26. A. Abulencia, et al. (CDF), Phys. Rev. Lett. 97, 062003 (2006) 
|arXiv:hep-ex/0606027v l] , ibid 97, 214003 (2006) | a rXiv:hep-ex/0609040y l] . 

27. V. Abazov et al. (DO), Phys. Rev. Lett. 97, 021802 (2006) |arXiv:hep-ex/0603029| v2] . 

28. H. G. Moser and A. Roussarie, Nucl. Instrum. Methods Phys. Res., Sect. A 384, 491 
(1997). 

29. E. Follana, C. T. H. Davies, G. P. Lepage, and J. Shigemitsu, Phys. Rev. Lett. 100, 
062002 (2008,) |arXiv:0706J726V2 [hep-lat]. 

30. B. L Eisenstein et al. (CLEO), Phys. Rev. D 78, 052003 r2008). 1a?Xrv:0806. 21121 ^3 
[hep-ex] . 

31. J. Rosner and S. Stone, "Decay Constants of Charged Pseudoscalar Mesons," review for 
the Particle Data Group [7]. larXiv:0802.1043k 3 [hep-ex]. For updated results see P.U.E. 
Onyisi et q/.(CLEO). larXiv:0901.lT47l [hep-ex], and J. P. Alexander et a/.(CLEO), 
rarXiv:0901.1216 [hep-ex]. 

32. N. Tantalo, "Lattice calculations for B and K mixing, " [ arXiv:hep-ph/0703241| yl] . 

33. M. Beneke, G. Buchalla, A. Lenz and U. Nierste, Phys. Lett. B 576, 173 (2003) 
ihep-ph/030734 4|. 



34. B. Aubert et al. (BABAR), arXiviOSOS.lTO [hep-ex]: K. F. Chen et al. (BELLE) Phys. 
Rev. Lett. 98, 031802 (2007) |arXiv:hep-ex/0608039v 4 [hep-ex]. 

35. P. C. Bhat, and H. B. Prosper, "Multivariate Methods in High Energy Physics: The 
Neural Network Revolution, World Scientific, ISBN-13: 978-9810243470 (2005); see also 



http : / / neur alnet s . web . cer n . ch/Neur alNets / nnwinhep . html[ 



36. R. Itoh et al. (BELLE), Phys. Rev. Lett. 95, 091601 (2005) [hep-ex/0504"030] ; B. Aubert 
et al. BABAR Phys. Rev. D 71, 032005 (2005) | hep-ex/0411016] . 

37. B. Aubert, et al. (BABAR), |ar2avi0807.4226v2 [hep-ex]; H. Ishino et al. (BELLE), Phys. 
Rev. Lett. 98, 211801 (2007) |hep-ex/0608035]. 



70 



38. L. Wolfenstein Phys. Lett. 13, 562 (1964). 

39. M. Gronau and D. London, Phys. Rev. Lett. 65, 3381 (1990). 

40. Y. Grossman and H. Quinn, Phys. Rev. D 58, 017504 (1998). See also J. Charles, Phys. 
Rev. D 59, 054007 (1999) |hep-ph/9806468|; M. Gronau, D. London, N. Sinha and R. 
Sinha, Phys. Lett. B 514, 315 (2001) [hep-ph/0105308] . 

41. I. Dunietz, H. R. Quinn, A. Snyder, W. Toki and H. J. Lipkin, Phys. Rev. D 43, 2193 
(1991). 

42. B. Aubert et al. (BABAR) Phys. Rev. D 76, 052007 (2007) arXiv:0705.2157l [hep-ex]. 

43. A. Somov et al. BELLE, Phys. Rev. Lett. 96, 171801 (2006) [ar Xiv:hep-ex/0601024| . 

44. A. F. Falk, Z. Ligeti, Y. Mr, and H. Quinn, Phys. Rev. D 69, 011502 (2004) 



hep-ph/0310242|. 



45. |http://ckmfitter.in2p3.fr7| 

46. A. E. Snyder and H. Quinn, Phys. Rev. D 48, 2139 (1993). 

47. B. Aubert et al. (BABAR), | arXiv:hep-ex/0703008> 2] . 

48. A. Kusaka et al. (BELLE), Phys. Rev. Lett. 98, 221602 (2007) | arXiv:hep-ex/0701015] . 

49. H. J. Lipkin, Y. Mr, H. R. Quinn, and A. E. Snyder, Phys. Rev. D 44, 1454 (1991); M. 
Gronau, Phys. Lett. B 265, 389 (1991). 

50. P. Robbe (LHCb), Nucl. Phys. Proc. Suppl. 170, 46 (2007). 

51. J. Charles, A. Hocker, H. Lacker, F.R. Le Diberder, S. T'Jampens, 
|arXiv:hep-ph/0607246v l]. 

52. M. Gronau and D. Wyler, Phys. Lett. B 265, 172 (1991); see also M. Gronau and D. 
London, Phys. Lett. B 253, 483 (1991). 

53. A. Bondar and T. Gerson, Phys. Rev. D 70, 091503 (2004) [ arXiv:hep-ph/0409281| . 

54. B. Aubert et al. BABAR, Phys. Rev. D 77 111102 (2008) arXiy:0802J052, [hep-ex]; K. 
Abe et al. (BELLE), Phys. Rev. D 73, 051106 (2006) | hep-ex/0601032] ; K. 

Gibson (CDF). [a?Xiv:0809.4809 : B. Aubert et al. BABAR. [aFXiv:0807.2408^y 2 [hep-ex]; 
K. Abe et al. (BELLE), Phys. Rev. D 73, 051106 (2006) [hep-ex/0601032] ; G. 
Marchiori (BABAR) larXiv:0810.0502l 

55. D. Atwood, I. Dunietz and A. Soni, Phys. Rev. Lett. 78, 3257 (1997) |hep-ph/9612433] . 



71 



56. Y. Horii et al. fBELLEl. [arXiv:0804.2063V l [hep-ex]; B. Aubert et al. (BABAR) Phys. 
Rev. D 72, 032004 (2005) [hep-ex/0504047] . 

57. R.H. Dalitz, Phys. Rev. 94, 1046 (1954). 

58. A. Giri, Y. Grossman, A. Soffer, and J. Zupan, Phys. Rev. D 68, 054018 (2003) 
|arXiv:hep^ph/0303187j . 

59. B. Aubert et al. (BABAR), Phys. Rev. D 78, 034023 (2008) [arXwi0804.2089 [hep-ex], 
they also include the mode KsK'^ K~ . 

60. K. Abe et al. fBELLE). 'i arXiv: 0803. 3375 1 [hep-ex]: 

61. A. Bondar and Poluektov, Eur. Phys. J. C47, 347 (2006); A. Bondar, A. Poluektov, 
iarXiv:0801.0840., [hep-ex]. 

62. P. Naik (CLEO), |irXiv:0810.3666' [hep-ex]. 

63. J. L. Rosner, Phys. Rev. D 42, 3732 (1990). 

64. A. S. Dighe, I. Dunietz, H. J. Lipkin and J. L. Rosner, Phys. Lett. B 369, 144 (1996) 
[hep-ph/9511363j . 

65. T. Kuhr fCDF). [a?XiV:0710.1789l [hep-ex]). 

66. A. S. Dighe, I. Dunietz and R. Fleischer, Eur. Phys. J. C 6, 647 (1999). 

67. T. Aaltonen et a/.(CDF), Phys. Rev. Lett. 100, 161802 (2008) arXiv:0712.2397Vl 
[hep-ex]; updated results in D. Tonelh et al. (CDF), arXiv:0810.3229v 2 [hep-ex]. 

68. V. M. Abazov et al. (DO). larXiv:0802.2255l yl [hep-ex]. 

69. J. Charles (CKM Fitter) 

| http ://www.slac. stanford.edu/xorg/ckmfitter/TalkBuffer/StatusOfTheCkmMatrixCapri2008.pdf 
see also O. Deschamps (CKM Fitter) [a?Xrv:0810.3139l yl [hep-ph]. 

70. J. Blouw (LHCb). [a^iv:0710.5124V 2 [hep-ex]. 

71. S. Stone and L. Zhang. [aFxrv:0812.2832 v3 [hep-ph]. 

72. B. Aubertet et al. (BABAR), Phys. Rev. D76, 031102 (20071 [aFXw:0704.0522| y2 [hep-ex]. 

73. B. Aubert et aZ. (BABAR), Phys. Rev. Lett. 93, 131801 (2004) [ arXiv:hep-ex/0407057] . 

74. B. Aubert aZ., (BABAR). laFXiv: 0807.4226 1 [hep-ex]. 

75. S.-W. Lin et a/.(BELLE), Nature 452, 332 (2008). 

72 



76. M. Peskin, Nature 452, 293 (2008). 

77. C. Jarlskog, Phys. Rev. D 35, 1685 (1987). 

78. M. Trodden, Rev. Mod. Phys. 71, 1463, (1999). 

79. R. Tarrach, Nucl. Phys. B 183, 384 (1981); A. S. Kronfeld, Phys. Rev. D 58, 051501 
(1998); P. Gambino and P. A. Grassi, Phys. Rev. D 62, 076002 (2000). 

80. I. I. Bigi, M. A. Shifman, N. G. Uraltsev and A. I. Vainshtein, Phys. Rev. D 50, 2234 
(1994); M. Beneke and V. M. Braun, Nucl. Phys. B 426, 301 (1994). 

81. N. Gray, D. J. Broadhurst, W. Grafe and K. Schilcher, Z. Phys. C 48, 673 (1990); 

K. Melnikov and T. v. Ritbergen, Phys. Lett. B 482, 99 (2000); K. G. Chetyrkin and 
M. Steinhauser, Nucl. Phys. B 573, 617 (2000). 

82. A. H. Hoang et al, Eur. Phys. J. C 3, 1 (2000). 

83. I. I. Bigi, M. A. Shifman and N. G. Uraltsev, Ann. Rev. Nuc. Part. Sci. 47, 591 (1997). 

84. P. Gambino, private communication (2006). 

85. A. H. Hoang, Z. Ligeti and A. V. Manohar, Phys. Rev. Lett. 82, 277 (1999). 
A. H. Hoang, Z. Ligeti and A. V. Manohar, Phys. Rev. D 59, 074017 (1999). 

86. E. Gardi for HFAG 2008 averages. 

87. CMcNehe, C. Michael and G. Thompson, Physics Letters B 600, 77 (2004). 

88. A. Renin and M. Steinhauser, Phys. Lett. B 538, 335 (2002). 

89. J.H Kuhn et al, Nucl. Phys. B 778, 192 (2007). 

90. A. H. Hoang, Phys. Lett. B 483, 94 (2000). 

91. N. Isgur and M. B. Wise, "Heavy Quark Symmetry," page 231, in B Decays Revised 
2nd Edition, ed. S. Stone, World Scientifc, Singapore (1994) . 

92. C. Glenn Boyd, B. Grinstein and R.F. Lebed, Phys. Lett. B 353, 306 (1995). 

93. J. Chay, H. Georgi and B. Grinstein, Phys. Lett. B 247, 399 (1990). 

94. C. Bernard et al. Submitted to Phys. Rev. D, larXiv:0808.2519J . 

95. N. E. Adam et al. (CLEO), Phys. Rev. D 67, 032001 (2003) [arXiv:hep-ex/0210040] . 

96. J. E. Duboscq et al. (CLEO), Phys. Rev. Lett. 76, 3898 (1996). 



73 



97. I. Adachi et al. (BELLE). [arXiv:0810.1657] [hep-ex]. 

98. B. Aubert et al. (BABAR), Phys. Rev. D 74, 092004 (2006). 

99. ICHEP08 updates available online at: 

|http://www.slac. stanford.edu/xorg/hfag/semi/ichep08/home.shtml 



100. D. Abbaneo et al. [ALEPH, CDF, DELPHI, L3, OPAL, SLD], |arXiv:hep-ex/0112028| See 
references therein. 

101. B. Aubert et al. (BABAR), Phys.Rev.D 77, 032002 (2008). 

102. B. Aubert et al. (BABAR), arXiv:0712.3493 ' [hep-ex] 



103. B. Aubert et al. (BABAR), arXiv:0809. 0828 [hep-ex] 

104. K.Abe, et al. (BELLE), Phys. Lett .B 526, 247 (2002). 

105. R. A. Briere et al. (CLEO), Phys. Rev. Lett. 89, 081803 (2002). 

106. I. Caprini, L. Lellouch and M. Neubert, Nucl. Phys. B 530, 153 (1998); C.G. Boyd, B. 
Grinstein, R.F. Lebed, Phys. Rev. D 56, 6895 (1997). 

107. B. Aubert et al. (BABAR), arXiv:0807.4978 ' [hep-ex]. 



108. B. Aubert et al. (BABAR), arXiv:0809.0 828 [hep-ex]. 

109. B. Aubert et al. (BABAR), arXiv:0712.3503l [hep-ex]. 

110. K. Abe, et al. (BELLE), Phys. Lett. B 526, 258 (2002). 

111. J. Bartelt, et al. (CLEO), Phys. Rev. Lett. 82, 3746 (1999). 

112. N. Isgur and M. Wise, Phys. Rev. D 93, 819 (1991). 

113. M. Okamoto (Fermilab Lattice, MILC and HPQCD Collaborations), Int. J. Mod. Phys. A 
20, 3469 (2005). 



114. B. Aubert et al. (BABAR), arXiv:0809. 0828 [hep-ex]. 

115. A. Anastassov et al. (CLEO), Phys. Rev. Lett. 80, 4127 (1998). 

116. V. Abazov et al. (DO), Phys. Rev. Lett. 95, 171803 (1995). 

117. B. Aubert et al. (BABAR), arXiv: 0708. 1738) [hep-ex]. 

118. D. Liventsev, T. Matsumoto et al. (BELLE), Phys. Rev. D 72, 051109(R) (2005). 



74 



119. K. Abe et al. (BELLE), Phys. Rev. Lett. 94, 221805 (2005); 
K. Abe et al. (BELLE), Phys. Rev. D 69, 112002 (2004). 

120. C. W. Bauer et al, Phys. Rev. D 70, 94017 (2004). 

121. O. Buchmuller, H. Flacher, Phys. Rev. D 73, 73008 (2006). 

122. K. Wilson, Phys. Rev. 179, 1499 (1969). 

123. D. Benson et al, Nucl. Phys. B 665, 367 (2003). 

124. M. Gremm and A. Kapustin, Phys. Rev. D55, 6924 (1997). 

125. A. Falk, M. Luke, and M.J. Savage, Phys. Rev. D 53, 2491 (1996); A. F. Falk and 
M. Luke, Phys. Rev. D 57, 424 (1998). 

126. P. Gambino and N. Uraltsev, Eur. Phys. J. C 34, 181 (2004). 

127. D. Benson et al, Nucl. Phys. B 710, 371 (2005). 

128. B. Aubert et al. (BABAR), Phys. Rev. D 69, 111103 (2004). 

129. C. Schwanda et al. (BELLE), Phys. Rev. D 75, 032005 (2007). 

130. A.H Mahmood et al. (CLEO) Phys. Rev. D 70, 32002 (2004). 

131. S. Chen et al, Phys. Rev. Lett. 87, 251807 (2001). 

132. J. Abdallah et al. (DELPHI), Eur. Phys. J. C 45, 35 (2006). 

133. D. Acosta, et al. (CDF), Phys.Rev. D 71, 051103 (2005). 

134. N. Isgur, Phys. Lett. B 448, 111 (1999). 



135. K. Abe et al. (BELLE), (ArXiv:hep-ex/0408139l ; ibid. [ArXiv:hep-ex/0409015| . 

136. A. H. Mahmood et al. (CLEO), Phys. Rev. Lett. 67, 072001, (2003). 

137. A. Mocker and V. KartvehshviU, Nucl. Instr. Meth. A 372 (1996) 469. 

138. P. Urquijo et al. (BELLE), Phys. Rev. D75, 032001 (2007). 

139. K. Abe et al. (BELLE), [arXiv:hep-ex/0611047| . 



140. K. Abe et al. (BELLE), [arXiv:hep-ex7() 508005 1 . 

141. P. Urquijo, "Measurement of the CKM elements \Vcb\ and \Vub\ from inclusive semilept 
B decays", PhD thesis (2008). 

75 



142. N. Uraltsev, Int. J. Mod. Phys. A 20, 2099 (2005). 

143. Presented by C. Schwanda, CKM08, Rome, Italy (2008). 

144. V. Aquila, P. Gambino, G. Ridolfi, and N. Uraltsev, Nucl. Phys. B 719, 77 (2005); A. 
Pak and A. Czarnecki, Phys. Rev. D 78, 114015 (2008); M. Bowling et al. Phys. Rev. D 
78, 074029 (2008). See also K. Melnikov. .arXiv:0803.0951v l [hep-ph]. 

145. K. Melnikov. laFxrv:0803.09'51l [heo-ph] (2008). 

146. T. Beicher et al. larXiv:0710.0680 i f2007l. 

147. B. M. Dassinger, T. Mannel and S. Turczyk, JHEP 0703, 087 (2007) 



arXiv:hep-ph/0611168|. 



148. N. E. Adam et al. (CLEO), Phys. Rev. Lett. 99, 041802 (2007) [ arXiv:hep-ex/0703041 1 . 



149. B. Aubert et al. (BABAR), Phys. Rev. Lett. 98, 091801 (2007) tarXiv:hep-ex/0612020] . 



150. T. Hokuue et al. (BELLE), Phys. Lett. B 648, 139 (2007) |arXiv:hep-ex/0604024| . 



151. B. Aubert et al. (BABAR), Phys. Rev. Lett. 97, 211801 (2006) tarXiv:hep-ex/0607089] 

152. K. Abe et al. (BELLE), [arXiv:hep-ex/0610054| . 

153. N. Isgur, D. Scora, B. Grinstein and M.B. Wise, Phys Rev. D 39, 799 (1989). 



154. P. Ball and R. Zwicky, Phys. Rev. D 71, 014015 (2005) |: arXiv:hep-ph/0406232] . 

155. D. Becirevic and A. B. Kaidalov, Phys. Lett. B 478, 417 (2000) |arXiv:hep-ph/9904490] . 

156. R. Hill, Phys. Rev. D 73, 014012 (2006) | arXiv:hep-ph/0 505129|y2] . 

157. K. C. Bowler et al. [UKQCD Collaboration], Phys. Lett. B 486, 111 (2000) 
| arXiv:hep-lat/9911011 | . 

158. J. Shigemitsu, S. Collins, C. T. H. Davies, J. Hein, R. R. Horgan and G. P. Lepage, Phys. 
Rev. D 66, 074506 (2002) | arXiv:hep-lat/0207011j . 

159. M. Okamoto et al, Nucl. Phys. Proc. Suppl. 140, 461 (2005) [arXiv:hep-lat/0409116] . 

160. E. Dalgic, A. Gray, M. Wingate, C. T. H. Davies, G. P. Lepage and J. Shigemitsu, Phys. 
Rev. D 73, 074502 (2006) [Erratum-ibid. D 75, 119906 (2007)] |arXiv:hep-lat/0601021] . 

161. M. C. Arnesen, B. Grinstein, I. Z. Rothstein and I. W. Stewart, Phys. Rev. Lett. 95, 
071802 (2005) faFXiv:hep-ph/0504209|. 



76 



162. J. M. Flynn and J. Nieves, PoS LAT2007, 352 (2007) larXiv:0711.3339] [hep-lat]. 

163. J. Bailey et al. ,arXiv:0811.3640l [heo-lat]. 



164. B. Aubert et al. (BABAR), Phys. Rev. Lett. 97, 211801 (2006) | arXiv:hep-ex/0607089] 

165. B. O. Lange, M. Neubert and G. Paz, Phys. Rev. D 72, 073006 (2005) 
arXTv:hep-ph/0504071| . 



166. J. E. Bartelt et al. (CLEO), Phys. Rev. Lett. 71, 4111 (1993). 

167. N. Isgur, J. Phys. G 18, 1665 (1992). 

168. 1. 1. Y. Bigi, M. A. Shifman, N. G. Uraltsev and A. L Vainshtein, Int. J. Mod. Phys. A 9, 
2467 (1994) (arXiv:hep-ph/9312359] . 

169. M. Neubert, Phys. Rev. D 49, 3392 (1994) [arXiv:hep-ph/9311325] . 

170. A. Bornheim et al. CLEO Collaboration Phys. Rev. Lett. 88, 231803 (2002). 

171. A. Limosani et al. (BELLE), Phys. Lett. B 621, 28 (2005). 

172. B. Aubert et al. (BABAR), Phys. Rev. D 73, 012006 (2006). 

173. I. L Y. Bigi, R. D. Dikeman and N. Uraltsev, Eur. Phys. J. C 4, 453 (1998) 
iarXiv:hep-ph/9706520j . 



174. C. W. Bauer, Z. Ligeti and M. Luke, arXiv:hep-ph/0111387 



175. C. W. Bauer, Z. Ligeti and M. E. Luke, Phys. Lett. B 479, 395 (2000) 
|arXiv:hep-ph/0 00216l] . 

176. B. Aubert et al. (BABAR), Phys. Rev. Lett. 100, 171802 (2008) arXiv:0708.3702l [hep-ex]. 

177. B.O. Lange, M. Neubert and G. Paz. Phys. Rev. D 72, 073006 (2005). 

178. J.R. Andersen and E. Gardi. JHEP 0601, 097 (2006). 

179. P. Gambino, P. Giordano, G. Ossola and N. Uraltsev, JHEP 0710, 058 (2007) 
.arXiv:0707.2493 [hep-ph]. 

180. V. Aquila, P. Gambino, G. Ridolfi and N. Uraltsev, Nucl. Phys. B 719, 77 (2005) 
jarXiv:hep-ph/0503083i . 

181. E. Gardi, JHEP 0404, 049 (2004) [arXiv:h ep-ph/04 03249l . 

182. U. Aglietti, G. Ferrera and G. Ricciardi, Nucl. Phys. B 768, 85 (2007) 
|arXiv:hep-ph/0608047l . 

77 



183. J. L. Rosner et al. (CLEO), Phys. Rev. Lett. 96, 121801 (2006) [arXiv:hep-ex/0601027] . 

184. M. B. Voloshin, Phys. Lett. B 515, 74 (2001) [arXiv:hep-ph/0106040] . 

185. N. E. Adam et al. (CLEO), Phys. Rev. Lett. 97, 251801 (2006) |arXiv:hep-e x/0604044] . 

186. J. Yelton et al. (CLEO), |arXiv:0903.060T] . 

187. M. Steinhauser. llarXiv:0809.T925 [hep-ph]. 

188. A. Gray et al. [HPQCD Collaboration], Phys. Rev. Lett. 95, 212001 (2005) 
|arXiv:hep-lat/0507015l . 

189. X. H. Guo and M. H. Weng, Eur. Phys. J. C 50, 63 (2007) |larXiv:hep-ph/0611301J. 

190. C. Bernard et al. [Fermilab Lattice, MILC and HPQCD Collaborations], PoS LAT2007, 
370 (2007). 

191. K. Ikado et al. (BELLE), Phys. Rev. Lett. 97, 251802 (2006). 

192. I. Adachi {et al.) (BELLE), arXiv:0809.3"8 34] . 

193. B. Aubert et al. (BABAR), Phys. Rev. D 77, 011107 (2008). 

194. et al. f BABAR). larXiv:0809.4027l [heo-ex]. 

195. S. Trine, lar2£iYi081Q.3633 [hep-ph]. 

196. Flavour in the era of the LHC CERN, 15-17 May, 2006, Yellow Reports 
(CERN-2007-004) and reference therein. 

197. G. Buchalla, A.J. Buras, M. E. Lautenbacher, Rev. Mod. Phys. 68, 1125 (1996). 

198. D.Atwood, M. Gronau, A. Soni, Phys. Rev. Lett. 79, 185 (1997); B. Grinstein, Y. 
Grossman, Z. Ligeti, D. Pirjol Phys. Rev. D 71, 011504 (2005). 

199. M. Misiak et al., Phys. Rev. Lett. 98, 022002 (2007). 

200. T. Becher, M. Neubert, Phys. Rev. Lett. 98, 022003 (2007). 

201. S. J. Lee, M. Neubert and G. Paz, Phys. Rev. D 75, 114005 (2007) 
|arXiv:hep-ph706 09224| . 

202. A. Limosani, in a private communication performed the average for this article. 

203. P. Gambino and M. Misiak, Nucl. Phys. B 611, 338 (2001) | arXiv:hep-ph/01 04034]. 

204. S. Chen et al. (CLEO), Phys. Rev. Lett. 87, 251807 (2001). 

78 



205. K. Abe et al. (BELLE), Phys. Lett. B 511, 151 (2001). 

206. B. Aubert et al. (BABAR), Phys. Rev. D 72, 052004 (2005). 

207. B. Aubert et al. (BABAR), Phys. Rev. Lett. 97, 171803 (2006). 

208. B. Aubert et al. (BABAR), Phys. Rev. D 77, 051103 (2008). 

209. A. Limosani et al. (BELLE), presented at Moriond E.W. (2008). 

210. T. Hurt, E. Lunghi, W. Porod, Nucl. Phys. B 704, 56 (2005). 

211. M. Gronau, Y. Grossman, D. Pkjol, A. Ryd, Phys. Lett. 88, 051802 (2002); M. Gronau 
and D. Pirjol Phys. Rev. D 66, 054008 (2002). 

212. B. Aubert et al, (BABAR), Phys. Rev. D 72, 051103 (2004); M. Nakao et al. (BELLE), 
Phys. Rev. D 69, 112001 (2004). 

213. B. Aubert et al. (BABAR), Phys. Rev. D 70, 112006 (2004). 

214. M. Nakao et al. (BELLE), Phys. Rev. D 69, 112001 (2004). 

215. M. Beneke et al, JHEP 0506, 071 (2005). 

216. B. Aubert et al. (BABAR), arXiv:0805. 4796 (submitted to PRL). 

217. S. Nishida et aZ.(BELLE), Phys. Rev. Lett. 93, 031803 (2004). 

218. T. E. Coan et al. (CLEO), Phys. Rev. Lett. 86, 5661 (2001). 

219. B. Aubert et al. (BABAR), Babar Collaboration Phys. Rev. D bf 77, 032007 (2008). 

220. T. Aaltonen et al (CDF), Phys. Rev. Lett. 100, 101802 (2008). 

221. V. Abazov et al . (DO), Phys. Rev. D 76, 092001 (2007). 

222. B. Aubert et al. (BABAR), Phys. Rev. D 77, 032007 (2008). 

223. M.C. Chang et al BELLE collaboration Phys. Rev. D 68, 111101 (2003). 

224. R. F. Harr, (CDF), "Rare Decays of B and D Hadrons at CDF," larXiv:0810.3444] 
[hep-ex] . 

225. B. Aubert et al. (BABAR), Phys. Rev. D 73, 092001 (2006). 

226. I. Adachi, et al. (Belle). [a?Xrv:0810.0335l [hep-ex]. 

227. A. Aaltonen et al. fCDFl. [a?Xrv:0804.3908l (2008). 

79 



228. M. Iwasaki et al. (BELLE, Phys. Rev. D 72, 092005 (2005). 

229. B. Aubert et al. (BABAR), Phys. Rev. Lett. 93, 081802 (2004). 

230. A. Ali, E. Lunghi, C. Greub, M. Walker, Phys. Lett. B 66, 034002 (2002). 

231. B. Aubert et al. (BABAR), Phys. Rev. D 73, 092001 (2006). 

232. B. Aubert et al. (BABAR), Phys. Rev. Lett. 93, 081802 (2004). 



233. A. Ishikawa et al, Phys. Rev. Lett. 96, 251801 (2006) |arXiv:hep-ex/0603018| . 

234. M. Pierini, "Update on the unitarity triangle (UTFit)," presented at 34th Int. Conf. on 
High Energy Physics (ICHEP 2008), Philadelphia, PA, July 2008. 

235. E. Blucher and W. B. Marciano, "Kd, Ks, THE CABIBBO ANGLE, AND CKM 
UNITARITY," in the Particle Data Group reviews [7]. 

236. This type of analysis was done with somewhat different assumptions earlier by K. 
Agashe, M. Papucci, G. Perez and D. Pirjol hep-ph /0509117] . 

237. Measurement of Asl for Bg mesons are given by V. Abazov et al. (DO), Phys. Rev. Lett. 
98, 151801 (2007); CDF Note 9015 07-08-16 available at 

[http://www-cdf.fnal.gov/physics/new/bottom/bottom.htmll See also S. Laplace, Z. 
Ligeti, Y. Nir, and G. Perez, Phys. Rev. D 65, 094040, (2002) [hep-ph/0 202"oTo] and 
references contained therein. 

238. There are considerable theoretical uncertainties, see A. Lenz and U. Nierste, JHEP 06, 72 



(2007) [hep-ph/0612167| 



80 



